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ABSTRACT 

Context. Hatchell et al. (2005) (Paper I) published a submillimetre continuum map of the Perseus molecular cloud, detecting the 
starless and protostellar cores within it. 

Aims. To determine the evolutionary stage of each submm core in Perseus, and investigate the lifetimes of these phases. 
Methods. We compile spectral energy distributions (SEDs) from 2MASS (1-2 /mi), Spitzer IRAC (3.6, 4.5, 5.8, 8.0 yt/m), Michelle 
(11 and 18 tars), IRAS (12, 25, 60, 100 /mi), SCUBA (450 and 850 jum) and Bolocam(1100 //m) data. Sources are classified star- 
less/protostellar on the basis of infrared and/or outflow detections and Class I/Class on the basis of 7b i, L\, \/L smm and F^^/Fgso. 
In order to investigate the dependence of these evolutionary indicators on mass, we construct radiative transfer models of Class 
sources. 

Results. Of the submm cores, 56/103 (54%) are confirmed protostars on the basis of infrared emission or molecular outflows. Of 
these, 22 are classified Class 1 on the basis of three evolutionary indicators, 34 are Class 0, and the remaining 47 are assumed starless. 
Perseus contains a much greater fraction of Class sources than either Taurus or Rho Oph. We derive estimates for the correlation 
between bolometric luminosity and envelope mass for Class I and Class sources. 

Conclusions. Comparing the protostellar with the T Tauri population, the lifetime of the protostellar phase in Perseus is 0.25-0.67 Myr 
(95% confidence limits). The relative lifetime of the Class and Class 1 phases are similar, confirming the results of Visser et al. (2002) 
in isolated cores. We find that for the same source geometry but different masses, evolutionary indicators such as 7b i vary their value. 
It is therefore not always appropriate to use a fixed threshold to separate Class and Class I sources. More modelling is required to 
determine the observational characteristics of the Class 0/Class I boundary over a range of masses. 

Key words. Submillimeter;Stars: formation;Stars: evolution;ISM: structure;ISM: dust, extinction 



1. Introduction 

The submillimetre ( submm) surveys of the Perseus molecular 
cloud wit h SCUBA (Hatchel feTaLlboOSL hereafter Paper I) and 



Bolocam ( lEnoch et al.ll200 6) have for the first time produced a 
fairly complete census of the star formation activity in this cloud, 
identifying all the submm cores above their flux detection limits. 
Previously, these submm cores have been discussed in terms of 
their masses, clustering and their relationship to the molecular 
cloud (lHatchell et al.ll200l lEnoch et alJlluOa iKirk et al.ll2006l) 
but so far there has been little attempt to identify their evolu- 
tionary status. With an identification of the protostars amongst 
the ~ 100 submm cores, questions of the star formation rate, 
the timescales of the prestellar and protostellar phases, and the 
detailed evolution of protostars can be addressed. 

The timescales involved in protostellar formation are partic- 
ularly important as constraints on the underlying physics con- 
trolling the protostellar collapse. If cores form at the point where 
turbulence rapidly completely dissipates, then collapse should 
occur on roughly a free-fall time. If, on the other hand, mag- 
netic fields and ambipolar diffusion play an important role in 
the evolution of the cor es, then the time for a core to collapse 
can be greatly increased (Ta ssis & Mou schovias 2004) . Current 
stimates for the ratio of the prestellar to the protostellar life- 



time vary from 2-3 dVisser et al.l2002l;lLee & Mversl 19991) to 1 
d Ward-Thompson et al.ll 1 994t I Jessop & Ward-Thompsonl2000l) . 

Estimates for the l ifetime of an embedded low-mass proto- 
star a r e 1-2 x 10 5 Myr dGreene et al.ll 1994tlKenvon & Hartmaniil 



lyr _ 

Il995t IWilking et alJI 19891 iKenvon et alJll990h . These early es- 
timates were based on infrared surveys (eg. IRAS) with very in- 
complete counts of protostars. Now that we have a submm-based 
census of the protostellar counts in Perseus, combined with 
much more complete estimat es of the T Tauri population down to 
below the brown dwarf limit dMuench et al.ll2003t IWilking et al.l 
l2004t |j0rgensen et alj 120061) . and new age estimates for the 
T Tauri population using updated pre-main-sequence (PMS ) 
evolutionary tracks ( IWilking et al. 200^ iMavne et al.llin prep]) , 
it is time to revisit these estimates. 

A catalogue of protostellar sources also provides a valu- 
able resource for followup study (eg. with interferometers 
such as ALMA) as subsamples with particular properties (star- 
less/protostellar, mass, luminosity) can be identified in order to 
investigate specific aspects of protostellar or starless core evolu- 
tion. 

An established method of differentiating between protostars 
and starless cores is on the basis of their spectral energy dis- 
tributions (SEDs). Perseus has now been surveyed not only in 
the mm/submm by SCUBA and Bolocam but also across the 
infrared by 2MASS, Spitzer and IRAS, so fluxes are available 
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across a wide range of wavelengths. Using these publicly avail- 
able data we can compile a set of spectral energy distributions 
for the submm sources in Perseus which is based on a homoge- 
nous dataset. These consistent spectral energy distributions are 
an important tool in determining evolutionary stage. 

In this paper we collate the spectral data on the submm 
sources in Perseus and discuss the evolution and lifetimes of pro- 
tostars within the cloud. In Sect. [2] we discuss the derivation of 
fluxes for the SEDs, including the new mid-infrared observations 
with Michelle on UKIRT, and in Sect. O the resulting SEDs. 
Sect[3]discusses the resulting classification and its implications. 
In Sect |4] we present models of early Class sources and their 
implications for the classification. In Sect.[3]we consider the life- 
time of the embedded protostellar phase and our conclusions are 
summarised in Sect. [6] 

Throughout, we assume a distance of 320 pc for the 
Perseus molecular cloud bas ed on the Hipparcos distance of 
IC34 8 (Ide Zeeuw et al.l 119991) . consistent wit h our earl i er pa- 
per (Paper J), tho ugh other recent studies dKirk et al.l 120061: 
lEnoch et"ail 120061) have assumed a closer distance of 250 pc 
based on extinction studies dCernis & Straizvsll2003l and refer- 
ences therein). 



2. Spectral energy distributions 

We have compiled a spectral energy distribution (SED) for each 
submm source in Perseus, based on an extended version of the 
850 /im catalogue given in Paper I. With a mass limit of of 
0.1 M Q for a 30 K protostellar envelope or 0.5 M Q for a 10 K 
starless core, the sources we expect to detect in the submm are 
those at the earliest stages of protostellar evolution - starless, 
Class and Class I. 



2.1. Contributions to the SED 

Spectral energy distributions were compiled from fluxes either 
from images or catalogues from the following datasets: 2MASS 
(1-2/im), Spitzer IRAC (3.6, 4.5, 5.8, 8.0 yum), Michelle (11 and 
18 /im) , IRAS (12, 25, 60, 100 u rn), SCUBA (450 and 850 fan, 
iPaper J) and Bolocam (1 100 //m, lEnoch et al.ll2006l) . 



2.1 .1 . SCUBA 850 /im and master catalogue 

SCUBA 850/mi flu xes are taken from our SCUBA map 
dHatchell et al.ll2005l) . The JCMT beam size at this wavelength 
is 14". 

The SCUBA chop results in poor sensitivity and noise arte- 
facts in reconstructed maps on angula r scales > 100" (see 
iJohnstone etaiteOOOt IvTsser et al.ll2002l Fig. 1). To avoid false 
detections and large uncertainties on measurements of the ex- 
tended flux we take out the large spatial scales by applying 
an unsharp mask (USM) filter. This filter removes structure on 
scales of > 2' (40, 000 AU at 320 pc). This was implemented 
by convolving the map with a 2' FWHM Gaussian and subtract- 
ing the resulting smoothed map from the original. This filtering 
also selects appropria te spatial scales for candidate protostars 
(eg. lMotte etai1ll998l) . 

To find integrated source flux es for the 850/im so urces 
we ran Cl umpfind (2D version) (IWilliams et al.l 1 19941) and 
Sextractor dBertin & Arnoutslfl996l) on the SCUBA data from 
Paper I. Clumpfind allocates pixels (and therefore flux) below 
the highest closed contour to sources using a friend-of-friends 
method. We also used Sextractor to extract fluxes within a fixed 



90" aperture. We ran Clumpfind with a base level of 3cr and in- 
crements of lcr (cr = 35 mJy/beam at 850//m and 200 mJy/beam 
at 450/im) and a minimum pixel count of 5 (3" pixels). 

The catalogue for which we calculate SEDs (the master 
catalogue) consists of (a) all sources with a peak flux above 
5cr at 850//m and a Clumpfind detection from the SCUBA 
850jum USM map; (b) Clumpfind detections with a peak flux 
above 3cr at 850/mi and either an identification in Paper I, 
or a Bolocam 1.1mm detection with Clumpfind (Enoch et 
al. 2005). We additionally included the well known Class I 
source IRAS 03410+3152 which is a 3<x 850/mi detection with 
Clumpfind in this category. The requirement for a Clumpfind 
detection means that we can allocate 850/mi source integrated 
fluxes to every source presented here in a consistent way. The 5<x 
detection limit provides a peak flux limited catalogue. Our aim 
in the wider compilation including the 3cr sources with multiple 
detections is to include as many real sources as possible while 
minimising the number of artefacts. 

The requirement of 5 pixels above 3<x with a 5cr peak should 
include any true point source with a peak above 5cr while ruling 
out noise spikes. The corresponding point source mass limit at 
850//m is 0.5 M at 10 K. The maximum mass within one beam 
of a starless core which escapes our detection limit is therefore 
0.5 M (see Sect. |2.2!2l for the details of the mass calculation). 
The masses calculated from flux within the Clumpfind area do 
not take into account low level flux < 3cr outside this area. 
Where the integrated flux is less than the peak flux per beam, 
the peak flux per beam is used to calculate masses. This affects 
12 low mass sources, mainly Bolocam detections with an 850//m 
peak < 5cr. The lower requirement of 5 pixels above 3<x with a 
peak of at least 3<x corresponds to a mass of 0.3 M Q at 10K. 

Several further sources which we identified by eye in lPaper j 
were also not identified by Clumpfind as separate sources but 
merged into one clump. Sextractor 90" aperture fluxes are also 
included in the SEDs for comparison with the Clumpfind fluxes 
to give an idea visually of the unc ertaintie s due to the differ- 
ing clump extraction methods. The iPapeTl source numbers (in 
some cases multiple) corresponding to each Clumpfind source 
are listed in Table [5] along with the cross-identification from the 
2MASS and Bolocam catalogues and other known names for 
each source. The source positions listed in Table [5] and used for 
cross-identification with sources at other wavelengths are the po- 
sitions of peak 850/mi flux. Pointing with SCUBA is typically 
good to within 3". 



2.1.2. SCUBA 450^m 

SCUBA 450jum integrate d fluxes (beam size 8") were also taken 
from our SCUBA maps dPaper j) . As for the 850jum data, 2D 
Clumpfind was run on the unsharp masked map. Contours were 
set to 3<t, 3<t + 2" x lcr with n = 0-6 and <x = 200mJy/beam. 
For many of the 850/zm sources, no corresponding 450//m peak 
was identified largely due to the higher relative noise level on 
the 450/mi map. The 450/mi fluxes are generally integrated over 
a smaller area than the 850//m fluxes due to the higher noise level 
at this shorter wavelength. 



2.1.3. Bolocam 



Bolocam 1.1mm fluxes were taken from lEnoch et alj d2006l) 
where Bolocam peak positions were within 20" of the 850//m 
peak. The large offset allowed takes account of the compara- 
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tively large beam of Bolocam (30" compared to SCUBA's 14"). 
The median offset between SCUBA and Bolocam peaks was 9". 

2.1.4. IRAS 

Upper limits on the IRAS fluxes were taken from HIRES 
processed maps. The HIRES processing was carried out with 
120 iterations at 60 and 100pm and 40 at 12 and 25pm, fol- 
lowing the method use d successfully on the L1448 cluster 
dWorf-Chase et al.l|2000|) . The resulting beams are typically 50- 
90" at 100pm and 25-60" at 12pm. Despite the HIRES process- 
ing many of the sources are still confused. As the 100pm beam 
sizes are usually large compared to the 850pm sources, the FIR 
flux for each submm core falls within one HIRES beam, and 
therefore the HIRES peak flux per beam for each source can be 
used as an estimate of the total flux. There may also be a contri- 
bution from nearby neighbours where sources are tightly packed, 
and for that reason the IRAS HIRES fluxes are marked as upper 
limits on the SEDs. 



2.1.5. Michelle on UKIRT 

Observations were carried out using MICHELLE on UKIRlQ 
on the nights of 14 September to 21 September 2002. The in- 
strument has a field of view of 67" by 50" with 0.21" pixels. 
Using 1 1.6pm and 18.5pm narrowband filters contiguous re- 
gions were imaged to cover the brighter SCUBA submillime- 
tre peaks, concentrating on the regions known to contain clus- 
ters: the IC348(HH211), L1448, L1455 and NGC1333 regions. 
During the observations the seeing had a FWHM of ~ 0.7" at 
1 1.6pm and ~ 1" at 18.5pm. Source fluxes were measured in 2" 
or 2.5" aperatures and the typical 1-sigma noise levels reached 
in the observations were 30 mJy at 1 1.6pm and between 90 mJy 
and 200 mJy at 18.5pm. 

2.1.6. Spitzer IRAC 

Spitzer IRAC maps at 3.6, 4.5, 5.8, 8.0 //m observed a s part 
of the Cores to Disks legacy project dEvans et al. 1 I2003I) were 
extracted from the archive using Leopard. Source positions 
and fl uxes were extracted using Sextractor dBertin & Arnoutsl 
I1996I) and source detection was carried out on the median map. 
Fluxes were measured within isophotal areas down to a thresh- 
old of twice the local background, which was typically 0.1- 
0.2 MJy Sr _1 at 3.6 pm and 6-15 MJy Sr -1 at 8 pm. Spitzer 
sources were considered to be associated with SCUBA sources 
where they lay within 12" or the Clumpfind source FWHM, 
whichever was the larger, of the peak. No presele ction was made 
on Spitzer colours, unlike J0r gensen et al.l d2006l) . 

2.1.7. 2MASS 

JHK fluxes were taken from the 2MASS point source catalogue 
for sources with rising (red) spectral indices H - K > 0.8 which 
lay within 12" or the submm source FWHM (as determined by 
Clumpfind), whichever was the larger, of the SCUBA 850pm 
peak. 



1 The United Kingdom Infrared Telescope (UKIRT) is operated by 
the Joint Astronomy Centre on behalf of the U.K. Particle Physics and 
Astronomy Research Council. 



2.2. Results 

SEDs for each source are shown in Fig. [7] Every SED includes 
the 850pm flux from Clumpfind, a Bolocam 1100pm flux or 
limit, IRAS HIRES upper limits, and Spitzer IRAC and 2MASS 
fluxes or upper limits. Additionally, an 850pm flux in a 90" 
aperture is shown for the majority of sources (those identified 
by SEXTRACTOR), a 450pm integrated flux for sources iden- 
tified by Clumpfind on the 450pm map, and Michelle 11.6 and 
18.5pm fluxes where available. For a few well-known sources 
further flux constraints are available in the literature but we have 
not included these as our aim was to compile as homogeneous a 
dataset as possible. The differing flux extraction methods, in par- 
ticular the necessarily different areas for flux measurement, re- 
sult in some mismatches in the SEDs. In a very few cases there 
are apparent mismatches between a submm source and an in- 
frared source with falling spectrum (sources 12, 13, 58, 62, 65, 
70). 

2.2.1. r bol andLboi 

For each source we calculated the bolometric temperature T\, \ 
and bolometric luminosity Lboi using a p iecewise integra t ion be - 
tween the flux measurements following iMvers & Laddl d!993l) . 
Where there is no Spitzer, Michelle or 2MASS detection, both 
7boi and Lboi are taken as upper limits to reflect the fact that the 
HIRES fluxes are upper limits. The bolometric temperature 7boi 
is robust to calibration uncertainties or missing fluxes with even 
a factor 2 change between submm and MIR wavebands chang- 
ing 7boi by only a few kelvin. The bolometric luminosities are set 
largely by the far-infrared limits set by IRAS HIRES, which we 
estimate in uncrowded regions are good to 30% but in crowded 
regions may lead to overestimates of the luminosity due to flux 
contributions from neighbouring sources. 

2.2.2. Masses 

We calculate envelope masses based on a constant dust temper- 
ature of 7d = 10 K and a 850pm dust opacity of 0.012 cm 2 g -1 . 
The temperature of 10 K was required by the constraints of 
the SCUBA 850pm a nd IRAS HIRES lOO^tn fluxe s on model 
greybodies (assuming Oss enkopf & Hennind d!994h icy coagu- 
lated dust type 5). Even at 10 K 13% of the 100 micron upper 
limits fall below the prediction of the greybody. Greybodies at 
r d = 10 K are plotted on the SEDs in Fig|7l note that the 850 mi- 
cron aperture fluxes (90" diameter aperture, solid triangles) are 
often lower than the Clumpfind fluxes and the 450 micron fluxes, 
which also often fall below the greybody, are calculated over a 
smaller area. Higher temperatures, even as high as 12 K, were 
ruled out for the main component of mass in most sources. For 
comparison, models of C lass and Class I sources dShirlev et al.l 
120021: lYoung et al.l l2003) suggest that the isothermal dust tem- 
perature with which to calculate representative masses should be 
13.8 +2.4 K for Class sources and 16+4 K for Class 1 sources, 
consistent with 10 K. Temperatures derive d from modelling the 
cores in Perseus as Bonnor-Ebert spheres dKirk et alj|2006l) in- 
dicate temperatures between 10-19K. The 10 K to which we are 
constrained by the greybodies lies at the low end of the expected 
temperature range, at least for the protostellar sources if not the 
starless cores. Thus our SEDs point to detection of a larger frac- 
tion of cold dust than is often measured for protostars, for which 
mean dust temperatures of up to 30 K are often assumed when 
calculating masses from submm fluxes. 
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A reason for the detection of more cold dust in our sample 
could be that the submm fluxes are taken from the Clumpfind de- 
convolution rather than more limited apertures, even though the 
map was prefiltered to remove structure on scales above 2' be- 
fore running Clumpfind. IRAS HIRES fluxes also provide tighter 
upper limits on the source fluxes than the IRAS PSC by reducing 
the contribution from nearby emission. However, the Clumpfind 
components are no bigger than the IRAS HIRES beam at lOOyum 
(typically 80") so the SED should give a fair representation 
of the clump. Considering all the uncertainties in calculating 
masses, we find that the true masses could be lower than our 
estimates but are unlikely to be higher. Firstly this is a result of 
the low Td compared to studies which assume higher tempera- 
tures - a factor of 3 compared to 15 K or 6 compared to 30 K. 
Secondly, we assume 320 pc whereas 250 pc is also an arguable 
distance for Perseus dCernis & Straizvsll2003l) . which would re- 
duce masses by a factor 1.9. Note, though, that our masses 
are co nsistent with the Bolocam study of Perseus (lEnoch et al.l 
120061) . which also assumed 10 K, once their lower assumed dis- 
tance is tak en into account. Our dust opac ity of 0.012 cm 2 g -1 , 
taken from lOssenkopf & Henning d 19941) assuming a gas/dust 
ratio of 161 by mass, is also at the low end of the range. Some 
authors prefer 0.02 cm 2 g -1 at 850 yum, which also would reduce 
masses by a factor 1.7. In summary, with a different choice of 
factors and fluxes measured in a more limited aperture, masses 
could be up to a factor of 10 lower than these we give. 

The results for masses along with T\, Q \ and L^i are given in 
Table |U 

3. Classification 

The shape of the SED only depends on the dust distribution rel- 
ative to the heating source and this depends on many parame- 
ters - envelope mass, w hich determine s the overall optical dept h 
dlvezic & Elitzuri 119971) ; dumpin ess dlndebetouw et all 120061) : 
orientation of any outflow cavity dWhitnev et alj 2003 ): and lo- 
cation of heating source on the near or far side of the clump. 
However, there are well known general trends with evolutionary 
stage, as the overall enve lope mass reduces and the size of the 
outflo w cavity increases dWhitnev et al.l [2003t lArce & Sargentl 
120061) . Bearing the caveats in mind, we apply a number of evolu- 
tionary indicators to classify of our objects based on their SEDs. 
As we are largely interested in statistics, individual misclassifi- 
cations due to an edge-on or pole-on orientation are less impor- 
tant as to some extent they cancel out across the population as a 
whole. What is important is that we use a consistent set of mea- 
surements across our sample, and therefore can directly compare 
between sources. 

We assume our sources fall into the usual classification 
scheme of starless core, Class or Class I protostar. Starless 
cores are those that do not contain a hydrostatic core (though 
some starless cores are protostellar cores which will go on to 
form one); observationally they are cold and show no evidence 
for embedded infrared, radio or outflow sources. The boundary 
between the two protostellar classes, Class and I, is defined by 
the point in evolution when the mass of the envelope equals the 
stellar mass. However, we have no way of measuring the stel- 
lar mass and instead must rely on criteria from modelling and/or 
observation to separate Class and Class I sources based on evo- 
lutionary indicators derived from observables. 

The first two evolutionary indicators we use are the stan- 
dard ones: bolometric temperature L DO i and the ratio of submm to 
bolometric luminosity, Lboi/L smm . T\, Q i is plotted against envelope 
mass M em in Fig.Q] Lboi = 70 K is used to draw a line between 



Class and Class I embedded protostars dChen et alJI 19951) . As 
our Tboi are often upper limits where there are no near- or mid- 
infrared detections, we have plotted the identified protostars and 
the sources with upper limits separately (the SCUBA 850/iin de- 
tection limits are also shown). The highest mass sources are cold 
Class sources, which we would expect to evolve to less mas- 
sive and warmer objects. 

The ratio of submm to bolometric luminosity is plotted in 
Fig. [2] For the submm luminosity we estimate L13 from the 
850//m integrated flux, which is available for all our sources, 
measured in a bandwidth of 50 GHz and corrected to 1.3 mm 
using F1.3 = Fgso x (850/1300) 3 5 assuming that dust opacity 
is optically thin and follows a power law with index /3 = 1.5. 
The bound ary between Class and Class 1 sources was cal- 
culated bv lAndre etafl d 19931) to lie at Lboi/^u = 20,000, 
but in order to exclude well known Class I sources such as 
L1455 IRS2 we are required to lower it to Lboi/Li.3 = 3,000. 
The upper limits from IRAS introduce some uncertainty in cal- 
culating Lbob and it ma y be that the differi ng analysis of the 
IRAS fluxes here and in lAndre et"aT] d 19931) is responsible for 
some the difference. However,a reduction of a factor of 4 in 
Lboi/Li.3 for the Cl ass 0/Class I boundary was also required by 
IVisser et al.1 d2002l) using a similar dataset and it seems likely 
that the Andre et all d 19931) limit was simply set too high on the 
basis o f the Class I sources known at the time. The lAndre et al.l 
d 19931) limit was consistent with a simple model with a single 
accretion rate which may not apply in all cases: the limit should 
increase to account for increased accretion rates at higher mass 
with models suggesting that the boundary lies somewhere be- 
tween Lboi K M em 10 (alternatively, oc L smm 10 when M env is cal- 
culated with a fixed T& as here) and Lboi K M em 2 ' oc L smm 20 
(ISmithll2000l) . Note that as envelope masses and submm lumi- 
nosities are both proportional to the subm m flux, this indicator 
is equ ivalent to plotting Lboi vs. M env (eg. lAndre & Montmerlel 
dT994h . 

A third measurement linked to the shape of the SED is the 
ratio of Spitzer IRAC to submm emission. This has the advan- 
tage of being reliable even when the IRAS data are confused. 
In Figf2]we also plot the 3.6/850 fim flux ratio vs. the 850 fim 
flux together with rough ranges for this ratio based on the evolu- 
tionary models of IWhitnev et al.l d2003l) . as a third evolutionary 
indicator. 

Statistics (mean and standard deviation) for the three evo- 
lutionary indicators are given in Table Q] The mean value of 
log 10 Lboi corresponds to 78 K for the total population in Perseus 
and 91 K for the 56 confirmed protostars (excluding the sources 
with upper limits on Tboi; see Sect B.ll l. As a temperature of 
70 K is used to divide the Class and Class I populations 
(see Sect. 13. 2\ these low values indicate already that the pro- 
tostellar population contains a significant fraction of deeply 
embedded sources. The mean values of the other indicators 
(log 10 (L bo i/Li. 3 ), log 10 (L 3 . 6 /L 8 go)) also support this. 

3. 1 . Protostars and starless cores 

Near- or mid-infrared detections by Spitzer and/or 2MASS iden- 
tify 54/103 sources (52 + 7%) as definitely protostellar. An ad- 
ditional two sources are identified as protostellar on the ba- 
sis of molecular outfl ows by our 12 CO o utflow survey of a 
subset of this sample dHatchell et alj [2006. Paper III). These 
sources are bl-bN/S (2) and NGC1333 SK31 (47). These were 
not detected by Spitzer, in the case of bl-bN/S apparently be- 
cause the driving sources are still too deeply embedded, in the 
case of NGC1333 SK31 (46) because any associated infrared 
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Fig. 1. Evolutionary indicators vs. envelope mass: 7boivs.M e nv. Left: confirmed protostars only. Right: all sources. Upper horizontal 
line: Class O/Class I boundary at T\, Q \— 70 K dChen et al. I fl995t) . Lower horizontal line: 30 K. Solid vertical line: SCUBA detection 
limit. Dashed diagonal line: Early Class models (Sect. [4]). Class I sources (see Sect. I3.2t are marked by circles (magenta). 



Table 1. Mean and standard deviation for each of the evolution- 
ary indicators log 10 (7boi/O, log 10 (L b oi/£smm), log 10 (F^. 6 /F sso) 
for (top) confirmed protostars; (bottom) all sources, as plotted in 
Fig. [2] 



Confirmed protostars (56) 


Mean 


Standard deviation 


iog 10 (r M /x) 


1.96 (9 IK) 


0.33 


log 10 (L bo ,/(3000 x L smm J) 


-0.0256 


0.617 


l°gl ( F 3.6/ F 850) 


-2.92 (0.012) 


1.24 


All sources (103) 


Mean 


Standard deviation 


iog 10 (r bo i/A:) 


1.89 (78K) 


1.04 


log 10 (L bol /(3000 x L smm )) 


-0.0431 


0.691 


l°glo( F 3.6/ F 85o) 


-3.28 


1.07 



source is masked by the scattered light from the bright neigh- 
bour NGC1333 ASR1 14 (45). This brings the total of confirmed 
protostars in Perseus up to 56. A number of sources are first 
identified as protostellar by this survey, namely numbers 14 
and 15 (near IC348), 36 (SE of L1448), and 58, 60 and 61 (in 
NGC1333). 



The remaining 47 sources have no positive detection of a 
protostellar indicator (IR, outflow) and therefore could be star- 
less. Further infrared, radio or outflow observations are needed 
to detect any further protostars in this remaining sample, partic- 
ularly at the low end of the mass range where any IR emission 
is likely to be faint and outflow masses low. Whereas the two 
luminosity-based indicators rely on IRAS upper limits and are 
not useful for the classification of the non-detections, F3.6/FS50 
on the other hand clearly places all sources with IR limits in ei- 
ther the Class or starless category (Fig. [2] third panel). 

Table [37X1 gives the counts of starless and protostellar cores 
for each subregion of Perseus. The percentage of starless cores 
varies between regions from 17 to 68% compared to the over- 
all mean of 46%. To investigate whether this variation is signif- 
icant, we calculate the probability of these percentages occur- 
ing by chance assuming the sources are drawn at random from 
a population which is 46% starless. The resulting probabilities, 
calculated from the binomial distribution, are given in the last 
column of Table The small clusters (L1448, L1455, Bl) 
contain small numbers each and are consistent with the cloud 
average. NGC1333, which has a high fraction of protostars, and 
IC348 and the distributed population of isolated cores, which 
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Fig. 2. Further evolutio nary indicators vs . envelope mass. Top: confirmed protostars only. Bottom: a ll sources. Left: Lboi/ L srnm , 
converting Lsgo to Li 3 (lAndre et al.ll 19931) . Right: F^.e/F^o with classification based on the models of IWhitnev et al.l d2003l) . The 
boundaries between classes are shown at representative values but in fact there is overlap between classes because of the effects 
of orientation. Class I sources (on the basis of 2/3 indicators) are marked as circles (magenta) and luminosity/flux upper limits are 
given for sources with no Spitzer IRAC detection. 



contain larger fractions of apparently starless cores, all show sig- 
nificant variation from the cloud average. In each of these cases 
the probability of being consistent with a 46% starless fraction 
is less than 3%. This suggests that in individual regions, there is 
significant variation in the star formation rate on timescales of 
order the lifetime of a starless or protostellar core. Alternatively, 



it is possible that the starless population in NGC1333 is under- 
estimated as starless cores are hard to identify in richly clustered 
regions, and this may be biasing the cloud average to a low frac- 
tion of starless cores. Excluding NGC1333, the percentage of 
starless cores rises to 55% in the rest of the cloud, from which 
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Table 2. Counts of starless and protostellar cores by region. The 
last column gives the probability (from the Binomial distribu- 
tion) of the counts occuring by chance assuming the sources are 
drawn randomly from a population with a starless fraction of 
47/103. 



Source 


Starless 


Protostellar 


% starless 


Prob. 


NGC1333 


10 


26 


28% 


0.013 


IC348 1 


13 


6 


68% 


0.026 


L1448 


1 


4 


20% 


0.20 


L1455 


2 


5 


29% 


0.21 


Bl 


1 


5 


17% 


0.13 


Isolated 


20 


10 


67% 


0.016 


Total 


47 


56 


46% 





Including the filament to the S of IC348 which contains HH21 1. 



the IC348 and the isolated population are less significantly dif- 
ferent. 

Comparing the number of submm cores in IC348 and 
NGC1333 to their PMS populations of 348 and 213 stars, respec- 
tively (Sect|5J, shows that IC348 is three t imes as rich in submm 
cores as NGC1333. Spitzer observations dJ0rgensen et al.ll2006l) 
also show that the ratio of Class II+III to Class I sources is higher 
in IC348 than other parts of the cloud. The reason for this is that 
IC348 has been forming st ars for longer (> 3 Mvr. lLuhman et all 
120031: iMuench et alJl2003h and has had time to build up a larger 
fraction of Class II/III sources. IC348 has certainly been active 
for longer than NGC1333 and also, possibly, than the remainder 
of the cloud, though it is harder to be definite on this because the 
Spitzer observations do not rule out a more spatially extended 
distributed pop ulation of young PM S stars, as recently discov- 
ered in Taurus dSle snick et al.l 120061) . Spitzer area selection on 
the basis of high A v biases towards younger sources still associ- 
ated with molecular cloud material. 

We find 53% of submm cores and 57% of protostars in large 
clusters (IC348 and NGC1333). The fraction for protostars con- 
firms our pre-classification result, based on the sub mm core s 
alone, of at least 40% of star formation large clusters dPaper J) . 
These fractions are also consistent with the > 50% star formation 
in large clusters f ound by 2MASS for Perseus/Orion/Monoceros 
Carpentei 20001) and ~ 60% by Spitze r for O rion/Ophiuchus 



( Allen et al. in prep.l) . |j0rgensen et al.l d2006l) find for the 



youngest Spitzer-selected Class I sources in Perseus that a some- 
what smaller fraction of 39% lie in clusters. The reason for the 
discrepancy with our results is that not all the Spitzer sources 
associated with submm peaks satisfy the c2d colour criteria. 
|j0rgensen et al.l din prepj) find that positional association with 
MIPS sources identifies more sources in clusters as protostars 
than IRAC colours alone. Taken together the results are all con- 
sistent with 40% of current star formation in Perseus taking place 
in the distributed population of isolated stars and small (< 100 
member) clusters. 

From the definite detections, the ratio of protostars to star- 
less cores is at least 47/56 = 0.8. Previous estimates for the star- 
less/protostellar ratio are hi gher, r anging from 2-3 Visser et al] 
d2002l):lLee & Mversl (fl999l) to 1 dW ard-Thompson et al.ll994t 
iJessop & Ward-Thompsonll2000l) . This would appear to suggest 
that our survey is incomplete to prestellar cores, although it is 
difficult to believe that we underestimate the number of starless 
co res by the factor of 12 necessary t o reach the ratio measured 
bv lJessop & Ward-Thompsonld2000l) . 

Assuming that the star formation rate in Perseus is constant, 
we estimate that the average lifetime of starless cores above our 



detection limit is 0.8 times the lifetime of the embedded pro- 
tostars. If more protostars are later identified among the appar- 
ently starless population, the relative lifetime of detectable star- 
less cores will be reduced. On the other hand, prestellar cores 
begin their life below our detection limit so the full timescale for 
prestellar development will be longer than the time over which 
we could detect them in these observations. Starless cores with 
a flattened central density distribution which peaks below our 
mass per beam sensitivity limit could still contain considerable 
(distributed) mass and remain undetected by our survey. We con- 
clude that prestellar timescales are roughly equal to or longer 
than protostellar lifetimes. 

3.2. Class and Class I sources 

Based on the three evolutionary indicators in Fig. [2] (7b i, 
Lboi/ii.3' ^3.6/^85o)> 16/103 sources are certain Class I objects 
(see Table |4j. A further six are likely Class Is based on two 
out of three evolutionary indicators bringing the total number 
of Class Is to 22. More spectral/spatial information is required 
in order to determine the true nature of a further five sources 
(58,76,84,88, and Bolo62), which are classified Class I on T\, \ 
but are classified as Class on the other two indicators. For these 
sources, the temperature is high but the estimated luminosity is 
too low to be consistent with a Class I identification. We assume 
these are Class 0. The classifications based on each evolution- 
ary indicator, and the majority vote, are given in Tableland the 
spatial distribution of the sources is shown in Fig|6] 

Thus, of the 56 confirmed SCUBA-detected protostars, the 
number of Class Is is only 22 out of 56 (39%). The remain- 
ing 34 confirmed protostars must be Class 0. On the basis of 
these counts, SCUBA-detected Class 0s are more common than 
Class Is in Perseus. 

That raises the question of how many Class I sources we 
fail to detect with SCUBA because they fall below our detec- 
tion limit. We can only detect 10 K envelope masses of 0.3 M , 
and are complete above 0.5 M Q . Class I sources begin by def- 
inition with an envelope mass equal to the stellar mass and the 
envelope mass subsequently reduces. For Class Is forming stars 
below about 0.6 M , therefore, we would not expect to detect 
the envelope for more than half of the Class I lifetime. 

We can, however, estimate the number of Class Is which we 
fail to de tect with SCUBA as Sp itzer is more sensitive to Class I 
sources. |j0rgensen et al.l (120061) make an independent estimate 
of the number of Class I sources in Perseus from the Spitzer sur- 
vey, which covers a similar area as SCUBA, on the basis of their 
infrared spectral index in the Spitzer IRAC and 2MASS bands . 
As the infrared spectral index classification dGreene et al . 1994) 
does not include a separate category for Class sources (as at 
that time they were not detected i n the infrared) both Cla ss and 
Class I sources are included in the [j0rgensen et all(l2006l) Class I 
count of 54. Of their sources, 14 lie outside the SCUBA map 
boundary and 29 coincide with SCUBA-detected p rotostars. The 
remaining 1 1 sources in the |j0rgensen et al1d2006l) Class I list lie 
within the SCUBA map but were not detected in the submillime- 
tre. These could be Class I sources with low envelope masses 
below the SCUBA detection limit. Including these Spitzer de- 
tections, the total number of Class I sources in our mapped area 
could therefore be as high as 22 + 1 1 = 33, bringing the total 
number of protostars in Perseus up to 56 + 1 1 = 67. The 34 con- 
firmed Class sources therefore still make up more than half of 
the total number of protostars. This fraction may be even higher 
if more of the IR non-detections later turn out to be protostellar 
through identification of molecular outflows or radio sources, as 
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these must also be Class 0: the limits on the infrared emission 
rules out Class Is on the basis of Fi^/F^o (Figj2j>. 

The ratio of Class 0s to Class Is in Perseus is quite unlike 
Ophiuchus and Tau rus clouds where Class I sources outnum- 
ber C lass 0s by 1 0: 1 (I Andre & Montmerlel 19941 iMotte & Andrei 
|2001|) and exceeds even Serpens where at least 25% of s ources 
are Class (5 Class 0s ident ified by iHurt & B arson vl d 19961) 
compared to 19 Class Is from lKaas et al.l d2004l) ). Although it 
is possible that there are some further undetected Clas s I sources 
which were not identified by |j0rgensen et al.l (H006) it is hard 
to imagine (given the sensitivity of the Spitzer c2d survey) that 
these outnumber the detections by 10: 1. 

One explanation for the large numbers of Class sources in 
Perseus could be that Perseus has undergone a recent burst of 
star form ation - perhaps triggered by 40 Persei, a B O. 5 star in 
Per OB2 dWalawender et al.l 120051; iKirk et al.l 120061) - and this 
has produced a large number of protostars which are still in the 
Class phase. Perseus has been forming stars from at least 3- 
4 Myr ago IC348 u ntil the present day dLuhman et al.l 120031: 
iMuench et al.l 120031) . Therefore multiple triggering events 
would be required to explain both the earlier formation of IC 348 
and the star formation in the last few xO.l Myr, which would be 
necessary to explain a current population of Class 0s. The aver- 
age star formation rate in Perseus is consistent with steady star 
formation and a protostellar core lifetime of a few xlO years, 
in agreement with previous studies (see Sect. [5] Paper I). 

An obvious difference between Ophiuchus and Taurus on the 
one hand and Perseus, Serpens and Orion on the other is that the 
protostellar envelopes are more massive in Serpens, Perseus and 
Orion. Roughly similar flux sensitivity limits mean observations 
trace the higher-mass end of the population in these more distant 
clouds (Perseus 320 pc, Serpens 259 pc, Orion 400 pc compared 
to Ophiuchus and Taurus, both at ~ 150 pc). Th e Class sources 
in Ser pens have envelope masses above 0.6 M Q dHurt & Barsonvl 
119961) : the Perseus Class 0s all have masses above 0.5 M (see 
Fig®. In fact, no Class sources with envelope masses less than 
~ 0.2 M are found in any of these clouds, o r indeed in the re - 
cent compilation of known Class so urces dFroebrichll2005l) . 
This is very clearly seen in Fig. 2 of iBontemps et all d!996l) . 
and is partly due selection effects: targetted submm surveys were 
largely based on known infrared and radio sources and therefore 
only picked up Class Is and a limited number of high mass (lumi- 
nous) Class 0s. Our SCUBA survey, on the other hand, is more 
sensitive to Class 0s than starless cores or Class I sources, though 
we may be missing the equivalents of many of the lower-mass 
objects in Rho Oph and Taurus in Perseus. Note that the majority 
of Perseus Class sources would still have masses above 0.2 M 
even if the masses were reduced by the maximum uncertainty 
factor of 10. 

The higher mass of the Perseus population could be directly 
responsible for longer Class lifetimes, if the evolution of proto- 
stars depends on mass in a way which favours long Class life- 
times for higher mass objects. Such a varia tion is predicted, for 
exam ple, by gravoturbulent fragmentation dSchmeja & Klessenl 
|2004|) . It is also a consequence of the competitive accretion 
mode l dBonnell et all 1200 U iBate et all 120031: iBonnell & Batel 
120061) . as higher mass cores continue to acquire material from 
the surrounding cloud during the Class phase and only em- 
bark on the Class I loss of envelope mass when the reservoir of 
material has run out. 

The Class I/Class classification of protostars is also af- 
fected by mass, because in practice the classification relies on 
observational indicators such as 7boi rather than a direct com- 
parison of the mass in the envelope and the stellar mass, which 
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Fig. 3. Histogram of masses for all sources (solid line), Class Is 
(dot-dash), Class 0s (dashed) and starless cores (dotted line). 



cannot be directly measured. More massive sources also have 
greater optical depth in their envelopes and therefore SEDs 
which peak at longer wavelengths for the same dust geometry. 
As a source evolves, loses envelope mass and evacuates a cav- 
ity, it emits more infrared radiation relative to submillimetre, and 
passes from Class to Class I on the basis of L smm /Lboi or ?boi- 
More massive sources, with greater initial optical depth, have to 
lose relatively more of their envelope mass to reach the required 
Ibd threshold to be classified as Class I. The lower envelope 
mass sources have higher 7b i and L smm /Lboi because of lower 
envelope optical depth and apparently make the transition from 
Class to Class I earlier in their evolution. Fig. [3] which shows 
the mass distribution for the cores separated by class, shows that 
the Class sources tend to have higher masses than the Class Is. 
In order to investigate this effect further, we model the emis- 
sion from a Class geometry envelope with a range of masses 
in Sect. |4] and discuss the implications for the classification of 
Class I and Class sources, but before we can do this we con- 
sider the relationship between mass and luminosity. 

3.3. Relationship between L bol andM em . 

The relationship between luminosity and envelope mass is 
shown in Fig. [4] Mass is measured from the submm flux assum- 
ing a constant temperature of 10 K whereas luminosity is mea- 
sured by integration across the entire spectrum and thus is most 
dependent on the IRAS fluxes. Note that although a higher lu- 
minosity source will have higher temperature components than 
a lower luminosity source of equivalent mass, the bulk of the 
mass is in cold gas therefore the constant temperature assump- 
tion in the mass calculation still holds. There is a general trend 
for higher mass sources to have higher luminosity. But there is 
not a tight relationship between envelope mass and luminosity. 
This scatter is intrinsic, and due to the fact that we are plot- 
ting not just Class I protostars (as iReipurth et al.l (Q993)) but 
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also Class Os, which have higher masses for th e same luminos- 
ity. Th is scatter was also seen in the sample of iBontemps et alj 
d 19961) . and it is this scatter which enabled us to use Lt, \/L smm 
which is proportional to Lboi/M env to estimate the evolutionary 
state of a core (Sect. [3). 

Higher mass accretion rates - and therefore higher lumi- 
nosities - are clearly driven by higher initial envelope masses. 
The evolution of luminosity with mass depends on the de- 
tails of the accretion process, and evolutionary tracks for Lb i- 
M env have been s uggested based on various accretion models (see 
lFroebrichll2005l and references therein). A possible scenario is 
that mass accretion rates and therefore accretion luminosity ini- 
tially increase exponentially and then reduce, resulting in M env - 
Lboi tracks which rise rapidly in luminosity (up and slightly to 
the left in Fig. [U during the Class stage and then decrease in 
both luminosity and mass (down and to the left in Fig.|4| >. An ex- 
ample of such an evolutionary track, from I Smith! d2000l Fig. 6), 
scaled by a factor 30 in mass in order to match the high envelope 
masses in Perseus, is plotted in Fig. |4] In this example the lumi- 
nosity decrease during the Class I phase follows Lboi K M em a 
with power law index a between 1 and 2. Of course, accretion 
rates and luminosity may also vary during both the Class and 
Class I evolutionary stages if accretion is variable, as is known to 
be the case for mass ejection in the jet (eg. lHartigan et al.l2005l) ). 
so evolutionary tracks are probably not as simple as this model 
represents. 

Luminosity evolution during the Class stage can produce 
a selection bias towards high mass Class 0s as they reach lu- 
minosities at which infrared or outflow emission can be de- 
tected earlier in their lifetimes. Therefore, higher mass Class 0s 
may be visible as protostars for longer. This bias is strong 
if the luminosity evolution is slow and weak if it is fast. If 



Class phase dHenriksen et al. 


119971 ISchmeia & Klessenll2004l 


IWhitworth & Ward-Thompson 


200 It IStamatellos et alj 120051). 



then the Class 0s will spend most of their lifetime close to the 
Class I boundary and the number of Class 0s which remain un- 
detected will be small. There is supporting evidence for this 
fast luminosity evolution in Fig. [4] as high mass sources with 
low luminosities are absent. This suggests either that the time 
taken to reach high luminosities (roughly, within two orders of 
magnitude of the maximum for a given mass) is much shorter 
than the subsequent evolution through the Class and Class I 
phases, or the proposed evolutionary track is wrong. For exam- 
ple, high mass Class sources could evolve from lower-mass 
Class sources through competitive accretion. 

We fitted a straight line to log 10 (Lb o i) vs. log 10 (M env ), as- 
suming uncertainties in log 10 (M env ) of 0.5 (factor of 3) and 
log 10 (Lboi) of 0.3 (factor of 2). For the entire sample of known 
protostars, 

log 10 (L bo i) = log 10 (M env ) x (1.96 + 0.36) - (1.09 ± 0.30). 

Fitting just Class sources, the gradient is similar but the lumi- 
nosities are lower: 

log 10 (L bo i) = log 10 (M env ) x (1.88 + 0.46) - (1.43 ± 0.46). 

Fitting only Class Is, the fit is 

log 10 (L b oi) = log 10 (M env ) x (1.85 + 0.52) - (0.46 ± 0.35), 

with a gradient consistent with the index of 1.6 found for 
Class I sources driving Herbig-Haro objects dReipurth et al.l 
119931) . These fits are plotted in Fig. [4] There are two selection 
effects: firstly, low-luminosity low-mass Class 0s are selected 
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Fig. 4. Bolometric luminosity vs. envelope mass for confirmed 
protostars. Class I sources (see Sect. [3]) are marked with circles. 
Power law fits to all sources, Class 0s and Class Is are shown 
(see text). An evolutionary track from ISmithl d2000l) . scaled by 
a factor 30 in mass, is shown for illustration with arrows at 
log 10 (age) = 3.5,4.0,4.5,5.0,5.5 yr (in green) 



against so the slope of the best fit line for Class 0s and the whole 
sample are probably too low; secondly, we believe some of the 
low mass Class I sources should be reclassified as Class 0s (see 
Sect. |4]i, so the slope of the best fit line for Class Is is proba- 
bly too high. The uncertainties on the fits are high and the scat- 
ter is large, so they should in any case be treated with caution. 
Nonetheless, the results are consistent with Class Is having gen- 
erally higher luminosities than Clas s 0s and a po wer law index 
between 1 and 2, consistent with the lSmithl (120001) model. 



4. Radiative transfer modelling of Class geometry 
sources with different masses 

We computed radiative-transfer models of Class obj ects us- 
ing the Monte-Carlo ra diative-transfer code torus dHarriesI 
l2000t lHarries et al.ll2.Q04b . The code i s based on the radiative- 
equilibrium method described bv lLucvl d 19991) and uses an adap- 
tive mesh in order to adequately resolv e sub-solar-radius scal es 
in parsec-scale computational domains dKurosawa et al.ll2004l) . 

Our Class geome try is based on the models described 
by IWhitnev et al.l d2003l) . In brief, this consists of a blackbody 
source of radiation surrounded by a flared accretion disc which is 
in turn embedded in an infalling envelope. Low-density bipolar 
cavities exist perpendicular to the disc. For simplicity we adopt 
ISM-type grains, although more sophisticated models could al- 
low for evolution of both the grain size distribution and the dust 
chemisty. We give full details of the density description, along 
with the dust grain size distribution and chemistry adopted, in 
the Appendix. 

4.1. Model parameters 

We computed models for a range of envelope masses between 
0.05 M and 10M o . The luminosity for each model was selected 
by adjusting R t to follow log 10 (Lboi) oc 1.5 * log 10 (A^env) with 
luminosi ties appropriate to mas ses calculated at 30 K (consis- 
tent with IBontemps et all d!996l) ). The resulting absolute lumi- 
nosities are higher (for the same mass) than for the Class 
sources (Sect. [3~3l ) but the relative scaling (and therefore the rel- 
ative effect on 7b i) is consistent. A constant M env :M» ratio of 3 
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M cnv /M G 


0.05 


0.1 


0.2 


0.5 


1.0 


3.0 


10.0 




67.4 


68.1 


68.1 


54.8 


48.1 


40.1 


34.9 



Table 3. Dependence of T\, \ with M env for Class models with 
identical geometry. 



dWhitney et al.ll2003h was assumed to represent the same evolu- 
tionary stage for sources of different masses. The basic model 
parameters are listed in Table lA.il 

We computed all the SEDs assuming an inclination of 90°; 
the shape of the SED is essentially independent of viewing angle 
except for pole-on orientations, in which case the near-to-mid 
IR portion of the spectrum is enhanced by direct and scattered 
radiation that escapes along the polar cavities. 

4.2. Model results 

The dependence of Tboi with M env for Class models with iden- 
tical geometry is given in Table 14.21 

This dependence is plotted as a dashed line in Fig[T] Low 
mass Class sources have significantly higher 7b i than their 
higher-mass counterparts, for the same geometry. Our models of 
young Class sources reach T^oi ~ 70 for 0.2 M Q and below and 
are on the borderline of being classified as Class I. Furthermore, 
although we chose to model a geometry representative of the 
very earliest Class sources, a 7b i-M env dependence also holds 
for other geometries as it is a fundamental property of radiative 
transfer, so models of more evolved envelopes with a higher lu- 
minosity to mass ratio and a wider outflow cavity would cross 
the Class 0/Class/I boundary at 7b i = 70 K at higher M env . Low 
mass sources are therefore more likely to be classified as Class I 
on the basis of T\, \ despite having a Class type envelope ge- 
ometry. 

4.3. Class and Class Is revisited. 

Using the models, we can revisit the classification of pro- 
tostars in Perseus and consider why we see so many more 
Class sources com pared to lower-mass samples such as 
IBontemps et al.l (Il996l) . All the Class sources which have 7boi 
above the predictions of the model as shown in Fig [TJ would be 
reclassified as Class I sources on the basis of 7boi if their en- 
velopes were reduced to < 0.2 M Q while retaining the Class ge- 
ometry. This means that essentially all our Class sources - ex- 
cluding only the three sources bl-bN/S (2), NGC1333 SK22(70) 
and HH340B(71) - would have been classified as Class I had 
they had the same geometry but lower mass. This conclusion 
holds even if we were to calculate envelope masses assuming a 
constant 30 K rather than 10 K, which reduces the mas ses by a 
factor of 7, for consistency with IBontemps et ail d 1996b : in that 
case 7 rather than 3 sources would retain their Class classifica- 
tion. 

Thus, if Perseus had been forming lower mass stars with en- 
velope masses under 0.2 M Q , but with the same envelope geom- 
etry, the fraction of sources classified Class t o Class 1 would 
drop to around 7/49, similar to the 20% found bv lBontemps et al] 
d 19961) . This potential misclassification of Class 0-geometry 
sources with low envelope massescould account entirely for the 
high number of sources classified Class in Perseus - which 
we believe are truly Class sources (see below) - compared to 
studies of lower mass regions where we believe many Class I- 
classified sources may in fact have Class geometry. 



If this is the case, and Perseus is representative, then the rel- 
ative lifetimes o f Class I and Class phases are not, as earlie r 
suggested, 10:1 dAndre & Montmerlelll994l; iMotte et alJI 19981) . 
but inst ead roughly 1:1, a s found in Perseus and in Lynds dark 
clouds dVisser et al.ll2002l) . assuming that the star formation rate 
is constant. 

How are the other SED-based evolutionary indicators shown 
in Fig |2] affected by mass? The SED-based indicator F^^/F^o 
should also decrease towards higher masses for the same en- 
velope geometry as the peak of the SED shifts longwards (see 
Fig. IA.2t . The remaining evolutionary indicator, L\, Q \/L smm , is 
based on input properties of the models as M env is derived 
from Ls mm . F rom Fig. |4j and as also suggested by models 
dSmithl l2000l) . the Lboi/M env boundary between Class I and 
Class sources should rise with M env and therefore so should 
Lboi/^smm- On this indicator, low mass sources should require 
lower Lboi/^smm (Lboi/M em ) to be classified as Class 0. The 
higher masses in our sample do not explain why we have to set a 
low L boi/Li 3 boundary at 3000 rather than 20,000 dAndre et al.l 
119931) . as we would expect Lboi/L smm at the Class 0/Class I 
boundary to rise with mass leading to a higher threshold for 
Perseus. 

Could many of the apparent Class sources in Perseus be 
misclassified and in fact be Class Is? As higher-mass sources 
have lower 7boi for the same geometry, it is the warm, high- 
mass sources in the Class classification which might in 
fact be Class Is despite Tboi less than 70 K. This assumes 
that the Class 0/Class I boundary lies roughly parallel to the 
M env -7boi track shown on Fig[T] only higher in Tboi- Therefore 
sources such as L1448 NW (27), NGC1333 IRAS2 (44), and 
NGC1333 IRAS7 (46), HH 211 (12), L1448 NW (27) and 
L1448 C (29) are the most likely sources to have a Class I-like 
geometry. In the case of some these sources their high Tboi can 
be explained by their multiplicity (see Sect J5.21 >. However, this 
list contains many classic bright Class sources, which cluster 
in this part of the M em -T\, \ diagram; if the number of Class 0s 
in Perseus is to be reduced by reclassification then Class pro- 
totypes such as HH211 and L1448 C will be among the first 
sources to go. This seems unlikely as the collimated molecular 
out flows from these sources also support the young age estimates 
(eg. iHirano et alj|2006t iDutrev et al.lll997l) . Therefore we con- 
clude that it is unlikely that a large number of our Class sources 
need reclassification on the basis of the mass dependence of the 
evolutionary indicators. It is more likely that we, too, have mis- 
classified low-mass sources with Tboi close to 70 K as Class I 
when they should be Class - eg. sources 84 and Bolo62 which 
fail on two evolutionary indicators. 

It is clear the evolutionary indicators - Tboi, L\, \/L smm 
F 3.6/^850 ~ which we use to classify protostars should all vary 
with mass. All of these evolutionary indicators are currently set, 
based on models or observations, at a fixed boundary appropri- 
ate to only a small range of masses. How exactly they vary at the 
Class 0/Class I boundary needs further investigation. It is clear 
that more radiative transfer modelling is needed in order to de- 
fine the observable characteristics of the Class 0/Class 1 bound- 
ary at M env = M* and before the numbers of sources in each class 
can be determined reliably. Alternatively, it may make sense to 
consider all the embedded sources (Class and I) as a single 
class with variation within that class. 

4.4. Protostellar brown dwarf candidates 

Our modelled Class geometry, with only a narrow outflow cav- 
ity, is intended to correspond to the earliest stage in protostellar 
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evolution. If this were true then these models would define a 
"birthline" for protostars with no protostars falling below the 
line. FigQ] shows that the majority of protostars do indeed lie 
above the line. Protostars must begin with lower luminosities 
than those we have modelled, but the early phase is both short- 
lived and difficult to detect in the infrared or outflows. Protostars 
then are predicted to e volve steeply upwards and to the left in 
the M env -r bo i diagram dSrruthll200dt lFroebrichll2005l eg. ). 

Protostellar brown dwarfs forming in isolation, if such things 
exist, will be found among the sources with envelope masses up 
to a few times the substellar mass limit and fitting the Class 
envelope model. Our lowest mass, lowest 7boi candidate Class 
brown dwarfs with Spitzer detections are sources 84 (0.77 M Q , 
74 K) and Bolo62 (0.5 M Q , 103 K), both of which are too mas- 
sive and warm to be serious Class brown dwarf candidates. 
Earlier, more embedded candidates may be found among the 
low mass sources without Spitzer detections, if any of these can 
be confirmed as protostellar. These sources include the Bolocam 
sources which were only detected by SCUBA at the 3<x level and 
cluster around the 3cr detection limit in FigQ] of which the low- 
est 7boi candidate is Bolo70, if this source can be demonstrated 
to be protostellar. Given the dependence of 7b i on M env and that 
in protostars the main source of luminosity is accretion, young 
brown dwarfs are likely to rapidly evolve to 7b i > 70 even in 
the Class phase. 

5. Lifetime of the protostellar phase 

To estimate the full duration of the protostellar phase, Class 
plus Class I, we can compare our census of protostars in Perseus 
to the older pre-main-sequence (PMS) T Tauri p opulation, for 
which age e stimates exist, follow ing the method of lKenvon et al.l 
d 19901) and IWilking et al.1 (Il989l) . Most of the Class II and III 
sources near the Perseus molecular clouds lie in the two main 
clusters IC348 and NGC1333. On larger scales, the OB asso- 
ciation Per OB2 contains 90% of the B s tars in Perseus but 
with photometric age estimates of 6—15 Myr dde Zeeuw & Brandl 
Il985t iGimenez & Clausen! 11994 we assume that IC348 and 
NGC1333 contain the only star formation in the last 5 Myr. 
Using the pre-main-sequence population and age estimates for 
these two clusters, we can estimate the recent star formation rate 
and compare to the currently forming population of embedded 
protostars. 

The most complete estimate of the population of IC348 is 
348 ± 47 stars with < 16.75 in an area of 20.5 x 20.5 arcmin 2 
dMuench et alj|2003l) . The mean age for t his cluster is 3-4 Myr 
dLuhman e t al. 2003} iMuench et al .1 1200*31 Mayne, priv.comm.). 
Any star formation before this time was either at a lower level 
or is an artefact due to the uncertainties in ages derived from 
colou r-magnitude diagrams dHartmannll200 U iBurningham et al.l 
\200± 

IWilking et al.l d2004l) detected a total of 213 stars in 
NGC1333 above a similar magnitude limit of niK = 16 (note the 
K-band luminosity function is falling off rapidly at these mag- 
nitudes, so the slight difference in sensitivity between the two 
PMS surveys makes little difference to the counts). In the centre 
of the northern cluster alone, they detect 1 26 stars of which an 
estimated 15 are field stars. Assuming the field star fraction is 
the same in the wider area, the total of pre-main-sequence stars 
is 188+25 stars. Ag e estimates for indiv idual stars in NG C1333 
vary f rom 0.3 Myr dWilking et al.ll2004l) to several Myr dAspinl 
120031) . though it is generally accepted to be less than 1 Myr old. 

The total number o f PMS stars in both c luste rs is therefore 
536 + 53 stars. As both lWilking et al.1 d2004l) and lMuenchet al.l 



d2003l) find that the IMF is falling off towards lower masses, bet- 
ter sensitivity to PMS stars within the observed cluster regions 
would not increase this number significantly; this is an almost 
complete census of the PMS population within these limited ar- 
eas. SCUBA, Spitzer and 2MASS observations all have found 
that the fraction of stars forming in the distributed population 
in Perseus, in isolated cores or small groups outside the main 
IC348 and NGC1333 clusters, comprises 20-40% o f the total 



(ie. an additional 20-60% over NGC1333 and IC348) dCarcenter 
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Hatchell et alj|2005 


lJ0rgensen et al.ll2006llMegeath & al. 


2006 


Allen et al. in prep. 


). To get the total past star formation 



rate in Perseus we correct the star counts in the clusters by this 
factor, including the ~ 20% uncertainties. An additional popula- 
tion of 20-60% increases the star count over the past 3.5 Myr to 
750+ 130 stars. 



We compare the PMS count to the population of 56 proto- 
stars identified in the submm (Sect. |3~2| i. 

In order to estimate a lifetime for the protostellar population 
we have to make an assumption about the current star formation 
rate. The simplest assumption is that the star formation has been 
constant over the lifetime of star formation in the region. There 
does not appear to be any strong evidence, from the Spitzer data 
or otherwise, that this is either a particularly active or inactive 
time for star formation in Perseus, which would justify using a 
different value than the time average for the current star forma- 
tion rate. Star formation in individual regions of the cloud does 
vary with time - for example, star formation in NGC1333 only 
began ~ 1 Myr ago. However, the temporal variations in the star 
formation rate in individual regions are assumed to average out 
when spatially averaged over the whole cloud. Clearly, the more 
regions one includes in such a calculation, the more likely this 
assumption is to hold, and in a single molecular cloud with just 
two main clusters and a distributed population, temporal varia- 
tions in the star formation rate are likely to be our greatest source 
of error. 

The average star formation rate derived from 750 ±130 PMS 
stars formed over the last 3.5 Myr is 214 stars per Myr with 95% 
confidence limits of 168-306 stars per Myr, where the uncertain- 
ties include only the uncertainties on the PMS count and age. 

A lower limit of 0.26 Myr on the protostellar lifetime is cal- 
culated by assuming the 56 confirmed submm-detected proto- 
stars are the entire protostellar population. However, this simple 
estimate does not take into account any corrections for either 
completeness and multiplicity, both of which may increase the 
true number of protostars above those detected. 

5.1. Completeness 

If we are going to compare source counts at the protostellar and 
PMS stages, we need to ensure that we are comparing like with 
like and selecting objects corresponding to a similar range of fi- 
nal stellar masses at each stage. The PMS studies extend well 
below the brown dwarf limit. Our SCUBA survey, on the other 
hand, is incomplete to low mass protostars: the 3<x detection 
limit for 10 K protostars is 0.3M Q . Spitzer c2d has a s imilar sen- 
sitivity limit to embedded sources dEvans et al.ll2003l) . though a 
complementary bias towards the more evolved cores, resulting 
in a small number of additional MIR detections of protostars. It 
is unlikely that either of these surveys would have detected the 
protostellar phase of any lower-mass PMS objects if they form 
individually. Stars of ~< 0.1 M may never have had envelope 
masses above the detection limit of 0.3 M Q and simply would 
not be detected by our SCUBA survey. As envelope mass re- 
duces with time until it reaches zero, essentially any protostar 
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ultimately drops below the submm detection limit. What frac- 
tion of the total number of protostars do we detect? The vari- 
ation of protostellar envelope mass with time is poorly under- 
stood, and made more complex by disparate estimates for the 
relative numbers of Class and Class I sources, so using models 
of envelope evolution does not seem to be a reliable basis for 
an estimate. The detection limits are also made more complex 
by the presence of protostellar disks, which can also be detected 
in the submm. Uncertainties in the envelope masses due to un- 
certain dust opacities, dust temperature and cloud distance in- 
troduce further systematic errors when comparing the two pop- 
ulations. In summary, we must take into account a fraction of 
undetected protostars, but exactly how large that fraction is is 
hard to estimate. 

5.2. Multiplicity 

We also have to correct for multiplicity. If a significant fraction 
of protostellar cores form multiple sources simultaneously, and 
those sources could be resolved in the infrared studies once they 
reached the PMS stage (ie. are not close binaries), then the num- 
ber of protostars is higher than the number of protostellar cores 
which we count and we will underestimate the protostellar life- 
time. Note that if cores form stars sequentially and then eject 
them then the number of cores will still represent the number of 
stars being formed at any one time, so sequential star formation 
does not introduce any error in our lifetime estimate. 

How much evidence is there for multiplicity in our sample? 
Some of the cores in our sample are known to be forming bi- 
naries or higher-order systems: could these result in multiple 
infrared detections? NGC1333 IRAS4A and B are each bina- 
ries or higher orderjystej^(|La^aD[i99l); NGC1333 SVS13 
is a binary dAnriada et"ai1 l2000l): L1448 IR S2 is a binary 
dWolf-Chase et al.1 l2000t lO'Linger et all 120061) . Most of these 
would not be resolved at the T Tauri stage in infrared studies 
equivalent to those of IC348 or NGC1333. The ov erall binary 
fraction in IC348 is 19 + 5% (iDuchene et al.ll 19991) and we es- 
timate the fraction of resolved bi naries to be < 20% of t his, as - 
suming the period distribution in iDuquennoy & Mayorl dl99ll) . 
On that basis 4%+l% of our sources will form binaries which 
would be resolved in the infrared at a later stage, which intro- 
duces only a small correction to the lifetime. 

Do we see any evidence for protostar-protostar multiplic- 
ity other than close systems which would not be resolved 
in the PMS studies? With a 14" beam we might expect to 
find confusion especially in the clustered regions, and we do. 
L1448 N A/B (27) is an obvious example. Bl-BN/S (2) is po- 
tentially a pair of protostars, though there is no clear detectio n 
of a second infrared source or outflow dHirano et al 1 ll999ah . 
Several NGC1333 sources are also multi ple when examined at 
higher resolution dSandell& Kneel 1200 ll) or with radio obser- 
vations iRodriguez et al.1 d 19971) : iRodrfguez & Reipurthl d!998l) : 
IRAS2A/B (44), SSV13 (43), IRAS7 (46), and ASR1 14 (45) are 
all believed to contain multiple protostars with separations that 
would be resolved in the infrared. These sources are among the 
most massive in our sample and therefore the easiest in which 
to detect multiplicity. On the above count of six sources we esti- 
mate 10% is a lower limit to the number of cores containing mul- 
tiple protostars, but the true number is likely to be more as only 
these massive, luminous protostars have been investigated with 
sufficient sensitivity and resolution. Note that this multiplicity is 
not easily identified in the infrared: many of these sources con- 
tain a luminous protostar which swamps fainter neighbours and 
in the few cases where one can resolve multiple Spitzer detec- 



tions, these appear to be due to shocked H2 outflow components 
rather than additional protostars. 

It is also possible that some of our cores contain both a proto- 
star and a collapsing prestellar core. Such sources would appear 
overmassive for their luminosity and therefore classified Class 0. 
Class sources account for slightly more than half our total pop- 
ulation of cores, which places an upper limit on the number of 
cores of this type. 

Increasing multiplicity also increases the protostellar life- 
time, as our submm core count again underestimates the true 
number of protostars. This effect is not entirely independent to 
that of incompleteness, however, as many of the objects formed 
in multiple cores will be low mass sources which would not have 
been detected as protostars if isolated because of their low enve- 
lope masses. 

5.3. Bayesian probability-based estimate of protostellar 
lifetimes. 

Our current state of knowledge does not allow us to place defi- 
nite numbers on either the mean fraction of the protostellar life- 
time to which our SCUBA survey is sensitive or the multiplicity 
of the submm sources. But it does allow us to place some limits 
and also make some informed estimates about the most likely 
scenarios. 

We can use this additional information to improve on the 
simple lower limit on the lifetime calculate earlier and estimate 
the range of possible protostellar lifetimes. To do this we use 
a Bayesian methodology to derive a probability distribution for 
the protostellar lifetime, based on probability distributions for 
the factors involved - the number of protostars, number of pre- 
main-sequence stars, and duration of star formation in IC348. 
The Bayesian treatment has the advantage of allowing arbitrary 
distributions to be treated rather than making the usual implicit 
assumption that all variables follow a normal distribution. In 
our case our state of knowledge about the number of proto- 
stars, which has a well-defined lower limit but a long tail to high 
numbers (allowing for high multiplicity or incompleteness), are 
poorly represented by a Gaussian, and therefore a more complex 
treatment is required. The other parameters - the age of IC348 
and the number of PMS stars - are represented by Gaussian vari- 
ables in the usual way. 

Considering in detail the probability distribution for the 
number of protostars, P(N proto ), we have a lower limit: we know 
that there are at least 56 protostellar cores. Including the addi- 
tional Spitzer detections, we take a count of 67 protostars to be 
more probable than 56. Taking into account the known multiplic- 
ity of at least 10%, we would anticipate that 74 protostars is even 
more probable. We represent this by a rising probability from 56 
to 74 following a partial Gaussian exp -(N ploto - 74) 2 /(2 x 9 2 ). 
Above this count we have to correct for the unknown number 
of undetected low-mass envelopes and the amount of multiplic- 
ity, about which our state of knowledge is very uncertain, but 
again we can place some limits. If protostellar envelope masses 
mirror the IMF, then we would expect our SCUBA detections 
above 0.3M H to acc ount for only 40% of the total population 
dMuench et al.ll2003l) and the total count could number as many 
as 140; however, as envelopes have, on average, masses greater 
than the resulting protostars and some lower-mass stars will form 
inside higher-mass multiple cores we consider such a high count 
to be unlikely. To represent this, we assign a probability above 
the peak probability at A^ ploto = 74 initially remaining high but 
ultimately falling steeply to a low but finite probability to val- 
ues close to the maximum of 140, following an exponential 
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Fig. 5. Probability distribution for protostellar lifetimes. Top: 
Assumed probability distribution for the number of proto- 
stars based on SCUBA detection of 56 above O.3M . Centre: 
Probability distribution for the star formation rate (SFR) based 
on the duration of star formation in IC348 and the number of pre- 
main-sequence stars in Perseus. Bottom: Resulting probability 
distribution for the protostellar lifetime. 



P(^V pl -o to ) °c exp -((A^proto - 74)/54.8) . The resulting compos- 
ite probability distribution is shown in the top panel of Fig. [5] 
These assignments of probabilities are necessarily somewhat ar- 
bitrary but consistent with the data and our state of belief about 
the system. As we will see, the large uncertainties in A^ ploto will 
manifest themselves as a broad spread in protostellar lifetimes. 

The probability distribution for the age of the clusters, P(A) 
we take as a simple Gaussian, mean 3.5 Myr, standard deviation 



0.5 Myr. The probability of the number of pre-main-sequence 
stars, P(A^pms), is also a Gaussian given by the product of the 
536 + 53 stars identified in NGC1333 and IC348 and a factor 
representing the 43% + 20% additional distributed population. 
The resulting distribution is also Gaussian with a mean of 652 
and a standard deviation of 97. From these two distributions - 
age and number of PMS stars - we derive the corresponding 
distribution for the star formation rate, shown in the middle panel 
of Fig. [5] using the formula: 



P(SFR) 



-f 

Jo 



P(A) P(JVpms - Ax SFR) dA. 



We then assume that the star formation rate has been constant 
since the onset of star formation in IC348 until the present day, 
and calculate the lifetime of the protostellar phase by comparing 
the star formation rate with the number of protostars, T proto = 
Nproto/SFR. The corresponding probability distribution, 
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is shown in the bottom panel of Fig. 

The result is that the most probable lifetime is 0.43 Myr with 
95% confidence limits of 0.25 - 0.67 Myr. The distribution is 
not Gaussian but skewed towards longer lifetimes, reflecting the 
tail of probability towards high numbers of protostars. There is 
a large range of lifetimes with high probabilities within the 95% 
confidence limits. The 95% confidence range interval is much 
more robust against changes in the details of the probability dis- 
tribution than the position of the peak probability. 

The uncertain submm detection rate for protostars and the 
uncertain multiplicity in the submm cores contribute the major- 
ity of the ~ 40% uncertainties in the protostellar lifetime. This 
will to a large extent be addressed by the next generation of 
more sensitive submm surveys, which have an order of magni- 
tude better mass sensitivity and will therefore detect envelopes 
down to below the substellar mass limit. Sensitive mid-infrared, 
outflow or radio observations will be required to determine if 
such low-mass objects are indeed protostars. The uncertain age 
spreads in IC348 and NGC1333 and the uncertain count of re- 
cent (< 4 Myr) T Tauri stars outside these main clusters also con- 
tribute to the uncertainties. Much recent work has been done on 
age sp reads in PMS clusters dHartmanrj200ltlBurningham et al.l 
120051) and these are unlikely to improve significantly but a better 
count of the distributed population may come from wide-field 
surveys. 

An embedded phase lifetime of 0.25-0.67 Myr (95% confi- 
dence limits) is consistent with the simple lower limit (0.25 Myr) 
we calculated at the beginning of this section but somewhat 
longer than the prev i ous estimates of 0.2 + 0.1 Myr in Rho Oph 
(IWilking et al.lll989HGreene et al.ll 1994 and 0. 15 + 0.05 Myr i n 



Taurus-Auriga dKenvon et al]l990tlKenvon & Hartmannl 1995b . 
Unlike these mid-infrared based studies, we are comparing with 
both a larger sample of T Tauri stars, due to more sensitive pho- 
tometry, and a much more complete sample of embedded stars, 
which includes Class objects to which the earlier surveys were 
blind, so it is unsurprising that we reach a different conclusion. 

6. Summary and conclusions 

The Perseus molecular cloud is a well known site of star for- 
mation containing a population of pre-main-sequence stars sim- 
ilar to Ophiuchus and intermediate between Taurus and Orion 
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dLuhman & Riekd 119991: iBriceflo et ail l2002t iLuhman et all 
120031) . We have compiled spectral energy distributions from the 
near-IR to mm wavelengths (using data from Bolocam, SCUBA, 
IRAS, Spitzer, Michelle (where available) and 2MASS) for the 
103 known submm cores across the molecular cloud. These 
cores lie within the boundaries of our SCUBA map of the 
Perseus molecular cloud, which covers roughly the area with 
13 CO in tegrated intensity greater than 4 K km s _1 or A v > 4 
dPaper J) . From the presence of a near- and/or mid-infrared 
source or molecular outflow, we have identified 56 cores which 
are definitely protostellar. A further 1 1 Class I sources have been 
identified by Spitzer, bringing the total number of definite proto- 
stars in Perseus to 67. There may be further protostars among 
our submm detections which are not yet luminous enough to 
be detecte d in the IR or be identified by their outflow activity 
dPaperffll) . 

We use three evolutionary indicators to classify the sources 
on the basis of their spectral energy distributions: 7b i, ^3.6/^850 
and L smm /Lboi. On the basis of these indicators, 22 protostars are 
Class I (33 including additional Spitzer-only detections) and 34 
are Class 0. This is a much higher fraction of Class sources 
than is found in Ophiuchus or Taurus. On the basis of these 
counts, the lifetimes for the Class phase are similar to the 
Class I phase, rather than much shorter as has been previously 
found. 

Protostellar envelopes in Perseus are more massive than 
those in Ophiuchus or Taurus and this is the main reason for 
the high number of Class 0-classified sources in Perseus. It is 
generally true that higher mass envelopes have higher optical 
depth envelopes and therefore SEDs which peak at longer wave- 
lengths, resulting in lower 7boi- We have demonstrated this by 
modelling a typical Class geometry at a range of different 
masses and find that the bolometric temperature rises as the mass 
decreases so that a source which would be classified Class 
at 1 M would be classified Class I below 0.2 M on the ba- 
sis of this evolutionary indicator. Low envelope mass sources, 
such as the majority in Taurus or Ophiuchus, are therefore more 
likely to be classified Class I whereas the higher mass envelopes 
which are detected in Perseus are more likely to be classified as 
Class 0. This classification issue arises because of the difficulty 
in defining the Class 0/Class I boundary in terms of observable 
evolutionary indicators, all of which (Tboi, L\, \/L smm Fi^/F&so) 
vary with envelope mass for the same geometry of source. More 
investigation with radiative transfer models is needed to deter- 
mine where the boundary should lie. However, despite the issues 
in classification we believe that the count of Class sources in 
Perseus is roughly correct and similar to the number of Class Is, 
and that therefore - assuming a constant star formation rate - 
the lifetimes of these two phases, at least for stars starting from 
envelopes of half a solar mass or more, are similar. 

There is also a possible physical explanation for why mas- 
sive clusters might contain more Class sources. If there is 
competitive accretion for material within the cloud, then more 
massive sources may continue to accrete material onto their en- 
velopes and maintain their Class status for longer. It is possible 
that this mechanism operates to maintain the Class sources in 
Perseus and that this is why there are so many in the clustered 
environments. 

We find that bolometric luminosity and envelope 
mass are roughly correlated with log 1() (Lboi/io) K 
log i n (M env / M ) 1 - 96±0 - 36 , consisten t with the results for Herbig- 



H aro flow s ( Reipurth et al.l fl9 93l) and the protostellar models 
of lSrnlthl (120001) . Class sources and Class I sources follow 



approximately the same power law but the Class sources have, 
on average, an order of magnitude lower luminosities. 

Our estimate for the lifetime of the entire embedded proto- 
stellar phase is 0.25-0.67 Myr. This estimate, higher than previ- 
ous estimates of ~ 0.2 Myr, is based on our survey of protostars 
with corrections for multiplicity and incompleteness, and with 
counts and age estimates for T Tauri stars from recent censuses 
of the populations in IC348 and NGC1333, assuming a constant 
star formation rate in Perseus over the last 3.5 + 0.5 Myr. 

The ratio of starless cores above the submm detection limit 
to protostars in Perseus is at most 0.8, though this could reduce 
as apparently starless cores may yet be identified as protostellar 
by more sensitive IR or outflow observations. The lifetime for a 
detectable starless core (> 0.3 M within one SCUBA beam) is 
therefore at least 0.8 times the protostellar lifetime or 0.4 Myr. 
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Fig. 6. Class (red c rosses), Class I (maroon stars) and starless cores (blue circles) overlaid on the Perseus molecular cloud outlined in C ls O (dark green contours 1 K km s => 
lHatchell et al.l d2005l) ) and 13 CO (light green contours 3 Kkm s _1 , Bell Labs 7m survey made available by John Bally http : //casa/colorado . edu/~bally). The area mapped 
by SCUBA is marked. 
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Fig. 7. Spectral energy distributions for the 103 submm detections in Perseus. Contributions from short to long wavelength: 2MASS 
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Fig. 7. Continued. 
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Table 4. Source properties: integrated 850yum flux Fgso; integrated 450/jm flux 
F 450 ; bolometric luminosity L bol ; envelope mass M cnv calculated at 10 K; the three 
evolutionary indicators - bolometric temperature r bo] ; L smra /Lbol; and F i6 /Fg S0 
(the number in brackets is the power of ten); and classification. Sources are 
classed as starless (S) if there is no infrared source or outflow, otherwise the 
classification is based on the three evolutionary indicators discussed in Sect[3] 
which, giving the condition for a Class I are: T\, \ > 70 K; Lb \/L imm > 3000; and 
^3.6/^850 > 0.003. The final classification is based on the majority of the indica- 
tors with results of individual indicators given in brackets in the order listed here. 
The commonly used source name is given in the last column; see[5]for alternative 
names and references. 



Source 


^850 


^450 


7*ol 


M clw 


^"bol ^bol/ ^smm 


F 3.6/^850 


Class 


Name 




Jy 


Jy 


L 


M e 


K 


x 1/3000 










1 


6.31 


37.8 


3.7 


17.7 


53 


0.5 


2.4( 


-5) 


(000) 


bl-c 


2 


9.31 


52.2 


< 2.5 


26 .1 


< 25 


< 0.2 


< 1.0( 


-5) 


(000) 


bl-bS 


4 


5.37 


50.4 


2.6 


15.0 


32 


0.4 


5.1( 


-5) 


(000) 


bl-d 


5 


3.67 


_ 


<0.3 


10.3 


< 33 


< 0.1 


< 2.6( 


-5) 


s 




7 


1.00 


_ 


1.3 


2.8 


158 


1.2 


1.4( 


-2) 


I (III) 


IRAS 03301+3057, Bl SMM6 


10 


0.32 


_ 


0.4 


0.9 


117 


1.1 


7.0( 


-3) 


I (III) 


Bl SMM11 


12 


8.19 


83.3 


4.5 


22.9 


31 


0.5 


l.8( 


-4) 


(000) 


HH211 


13 


3.59 


17.6 


1.9 


10.1 


59 


0.5 


4.6( 


-4) 


(000) 


IC348 MMS 


14 


1.98 


15.9 


5.6 


5.5 


180 


2.6 


5.8( 


-2) 


I (III) 




15 


1.57 


2.0 


1.9 


4.4 


67 


1.1 


9.1( 


-4) 


(010) 




16 


1.43 


7.2 


< 1.6 


4.0 


< 103 


< 1.0 


< 6.7( 


-5) 


s 




17 


1.18 


_ 


< 2.0 


3.3 


< 78 


< 1.6 


< 8.1( 


-5) 


s 




18 


0.61 


_ 


<0.6 


1.7 


< 86 


< 0.8 


< 1.6( 


-4) 


s 




19 


0.68 


_ 


< 1.2 


1.9 


< 61 


< 1.6 


< 1.4( 


-4) 


s 




20 


1.44 


_ 


<0.9 


4.0 


< 70 


< 0.5 


< 6.7( 


-5) 


s 




21 


0.62 


_ 


< 1.2 


1.7 


< 81 


< 1.8 


< 1.6( 


-4) 


s 




23 


0.32 


_ 


< 1.2 


0.9 


< 89 


< 3.3 


< 3.0( 


-4) 


s 




24 


0.56 


_ 


<0.9 


1.6 


< 93 


< 1.4 


< 1.7( 


-4) 


s 




25 


0.39 


_ 


< 1.0 


1.1 


< 85 


< 2.4 


< 2.5( 


-4) 


s 




26 


152.43 




<0.8 


2.7 - 


< 84 


< 0.7 


< 1.0( 


-4) 


s 




27 


11.03 




5.4 


30.9 


53 


0.5 


l.K 


-4) 


(000) 


L1448 NW 


28 


5.83 


196 


17.0 


16.3 


53 


2.7 


3.3( 


-4) 


(010) 


L1448 N A/B 


29 


6.51 


68.4 


6.3 


18.2 


49 


0.9 


l.K 


-3) 


(000) 


L1448C 


30 


5.85 


51.1 


4.1 


16.4 


50 


0.7 


l.2( 


-4) 


(000) 


L1448 IRS2 


31 


4.26 


26.5 


1.9 


11.9 


38 


0.4 


1.7( 


-5) 


(000) 




32 


4.36 


2.9 


<0.2 


12.2 


< 20 


< 0.0 


< 2.2( 


-5) 


s 




35 


1.89 


18.4 


7.0 


5.3 


105 


3.4 


5.3( 


-4) 


I (HO) 


L1455FIR4 


36 


1.91 


9.8 


3.1 


5.3 


59 


1.5 


5.6( 


-4) 


(010) 




37 


2.52 


20.5 


3.3 


7.1 


672 


1.2 


5.6( 


-2) 


I (III) 


L1455 PP9 


39 


0.72 


1.6 


0.8 


2.0 


394 


1.0 


7.3( 


-2) 


I (ioi) 


L1455 FIR1/2 


40 


0.89 


_ 


<0.6 


2.5 


< 63 


< 0.6 


< 1.1( 


-4) 


s 




41 


18.61 


131.0 


8.5 


52.1 


35 


0.4 


9.K 


-6) 


(000) 


NGC1333 IRAS 4A 


42 


9.00 


94.6 


7.6 


25.2 


39 


0.8 


8.9( 


-5) 


(000) 


NGC 1333 IRAS 4B 


43 


12.70 


159.9 


29.5 


35.6 


180 


2.1 


l.K 


-2) 


I (III) 


NGC1333 SVS13 


44 


9.73 


95.6 


22.7 


27.2 


58 


2.2 


2.9( 


-4) 


(010) 


NGC 1333 IRAS 2A 


45 


8.82 


91.3 


10.2 


24.7 


164 


1.1 


L.6( 


-2) 


I (III) 


NGC1333 ASR 114 


46 


6.10 


32.8 


6.1 


17.1 


67 


0.9 


6.4( 


-4) 


(000) 


NGC1333 ASR 32/33 


47 


4.33 


26.3 


< 1.8 


12.1 


< 54 


< 0.4 


< 2.2( 


-5) 


(000) 


NGC1333 SK31 


48 


2.92 


20.9 


3.2 


8.2 


49 


1.0 


2.0( 


-4) 


(000) 


NGC 1333 IRAS 4C 


49 


2.59 


21.1 


6.8 


7.3 


140 


2.4 


l.2( 


-2) 


I (III) 


NGC1333 SK6 


50 


2.58 


7.1 


5.6 


7.2 


88 


2.0 


2.7( 


-4) 


I (HO) 


NGC1333 HH 7-11 MMS 4 


51 


3.01 


10.1 


< 1.5 


8.4 


< 37 


< 0.5 


< 3.2( 


-5) 


s 


NGC1333 SK16 


52 


2.17 


_ 


4.4 


6.1 


57 


1.9 


3.5( 


-4) 


(010) 


NGC1333 HH7-11 MMS 6 


53 


2.14 


29.7 


< 3.7 


6.0 


<54 


< 1.6 


<4.5( 


-5) 


s 


NGC1333 SK26 


54 


3.01 


78.8 


106.9 


8.4 


138 


32.7 


1.6( 


-2) 


I (III) 


NGC1333 SK28 


55 


2.96 


14.4 


1.0 


8.3 


50 


0.3 


5.8( 


-6) 


(000) 




56 


1.96 




19.1 


5.5 


124 


9.0 


7.K 


-2) 


I (III) 


NGC 1333 SK29 


57 


2.10 


15.1 


< 8.0 


5.9 


< 81 


< 3.5 


< 4.6( 


-5) 


s 


NGC1333 SK33 


58 


1.56 




0.3 


4.4 


82 


0.2 


7.8( 


-5) 


(100) 




59 


2.21 




< 1.3 


6.2 


< 83 


< 0.5 


<4.3< 


-5) 


s 




60 


2.95 




< 1.2 


8.3 


< 215 


< 0.4 


< 3.3( 


-5) 


s 




61 


1.42 


2.5 


1.0 


4.0 


64 


0.7 


3.8( 


-4) 


(000) 




62 


1.21 


6.1 


2.7 


3.4 


50 


2.0 


7.K 


-5) 


(010) 


NGC1333 SK18 


63 


0.76 




2.8 


2.1 


103 


3.4 


5.K 


-2) 


I (III) 


NGC1333 SK32 


64 


0.92 




<0.3 


2.6 


< 52 


< 0.3 


< 1.0( 


-4) 


s 


NGC1333 Per4A3/4D 


65 


0.72 


2.1 


0.4 


2.0 


56 


0.5 


3.1( 


-3) 


(001) 


NGC1333 IRAS4B1 


66 


1.03 




<4.2 


2.9 


< 65 


< 3.8 


< 9.3( 


-5) 


s 


NGC1333 SK30 


67 


1.51 




8.4 


4.2 


127 


5.1 


1.8( 


-2) 


I (III) 




68 


1.34 


8.2 


1.3 


3.7 


69 


0.9 


7.4( 


-3) 


(001) 




69 


0.23 




0.04 


0.6 


436 


0.2 


l.8( 


-2) 


I (ioi) 




70 


0.53 




2.6 


1.5 


45 


4.5 


6.3( 


-4) 


(010) 


NGC1333 SK22 


71 


0.75 




0.03 


2.1 


37 


0.0 


2.5( 


-4) 


(000) 


NGC1333 HH340B 


72 


0.42 




<0.1 


1.2 


< 64 


< 0.3 


< 2.3( 


-4) 


s 
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Table 4. continued. 



Source 


^850 


F450 






Tb \ 


^bol/ ^smm 


F 3.6/^850 


Class 


Name 




Jy 


Jy 


L s 


M e 


K 


X 1/3000 








74 


0.24 


- 


0.06 


0.7 


144 


0.2 


5.4(-3) 


1 (ioi) 




75 


0.43 


- 


0.1 


1.2 


67 


0.3 


7.7(-5) 


(000) 




76 


5.60 


21.5 


2.4 


15.7 


109 


0.4 


4.2(-4) 


(100) 


Bl SMM5 


77 


2.26 


3.3 


0.7 


6.3 


60 


0.3 


7.7(-5) 


(000) 


Bl SMM7 


78 


1.49 


- 


4.7 


4.2 


181 


2.9 


2.3(-2) 


I (III) 


B5 IRS1 


79 


1.16 




< 1.4 


3.2 


< 175 


< 1.1 


< 8.3(-5) 


S 




80 


0.86 


12.1 


1.9 


2.4 


86 


2.0 


3.9(-3) 


I (III) 


IRAS 03235+3004 


81 


0.42 


- 


0.8 


1.2 


133 


1.8 


1.5(-2) 


I (HI) 


IRAS 03271+3013 


82 


0.65 




<0.7 


1.8 


< 53 


< 0.9 


< 1.5(-4) 


s 


Bl SMM1 


83 


0.19 




< 0.10 


0.5 


< 68 


< 0.6 


< 6.8(-4) 


s 




84 


0.28 




0.06 


0.8 


74 


0.2 


1.6(-4) 


(100) 




85 


0.11 




< 0.00 


0.3 


< 196 


< 0.1 


< 3.9(-3) 


s 


IRAS 03254+3050 


86 


0.79 




0.06 


2.2 


49 


0.1 


5.8(-5) 


(000) 




88 


0.55 




0.05 


1.5 


83 


0.1 


4.0(-5) 


(100) 




89 


0.66 


- 


<0.2 


1.9 


<45 


< 0.2 


< 1.4(-4) 


s 


Bl SMM9 


90 


0.33 




<0.9 


0.9 


<66 


< 2.7 


< 2.9(-4) 


s 


Per 7 


91 


0.58 


- 


<0.1 


1.6 


< 74 


< 0.2 


< 1.7(-4) 


s 


Per4A 


101 


0.13 




1.6 


0.4 


463 


35.2 


2.7(0) 


I (III) 


IRAS 03410+3152 


Bololl 


0.20 


- 


<0.2 


0.6 


<60 


< 1.0 


< 4.8(-4) 


s 




Bolo26 


0.13 


- 


< 0.01 


0.4 


< 64 


< 0.0 


< 8.7(-4) 


s 




Bolo27 


0.13 




< 0.05 


0.4 


< 103 


< 0.4 


< 9.4(-4) 


s 




Bolo44 


0.14 




< 1.0 


0.4 


< 55 


< 10.6 


< U(-3) 


s 




Bolo62 


0.18 


- 


0.07 


0.5 


103 


0.4 


4.0(-4) 


(100) 




Bolo64 


0.20 


- 


<0.1 


0.6 


< 237 


< 0.5 


< 4.8(-4) 


s 




Bolo70 


0.11 




<0.1 


0.3 


<42 


< 1.7 


< 1.5(-3) 


s 




Bolo89 


0.11 


- 


<0.9 


0.3 


< 58 


< 13.2 


< 1.5(-3) 


s 




Bolo90 


0.14 


- 


1.7 


0.4 


943 


24.9 


7.3(-l) 


I (III) 


IRAS 03380+3135 


Bolo92 


0.12 




< 1.1 


0.3 


< 68 


< 19.2 


< 1.9(-3) 


s 




Bolo94 


0.11 




< 1.1 


0.3 


< 71 


< 28.4 


< 2.8(-3) 


s 




Bololll 


0.18 




< 1.5 


0.5 


< 52 


< 7.6 


< 5.3(-4) 


s 




Bololl2 


0.15 




0.01 


0.4 


222 


0.1 


1.3(-2) 


1 (ioi) 




Bololl3 


0.13 


- 


<2.0 


0.4 


<59 


< 14.2 


< 7.4(-4) 


s 




Bololl4 


0.13 




< 0.01 


0.4 


< 56 


< 0.0 


< 7.2(-4) 


s 




Bolol21 


0.11 


0.9 


< 1.1 


0.3 


< 231 


< 34.2 


< 3.3(-3) 


s 




clf55 


0.23 




< 0.01 


0.6 


<67 


< 0.0 


< 7.8(-4) 


s 




clf68 


0.39 




< 0.03 


1.1 


<42 


< 0.1 


< 2.5(-4) 


s 




clf70 


0.51 


70.2 


<5.1 


1.4 


<34 


< 9.3 


< 1.9(-4) 


s 




clt74 


0.83 




<0.2 


2.3 


< 84 


< 0.2 


< 1.2(-4) 


s 




clf77 


0.50 


24.5 


< 1.1 


1.4 


< 26 


< 2.0 


< 1.9(-4) 


s 




clt79 


0.58 


3.7 


<0.2 


1.6 


< 27 


< 0.3 


< 1.7(-4) 


s 
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Table 5. Source peak positions a nd cross-identification with [Hatchell et alj 
j2005T) . Bolocam dEnoch et alj|2006l) . 2MASS, and other names commonly as- 
sociated with the source. 



Source 


RA(J2000) 


Dec(J2000) 


Paper I 


Bolocam 


2MASS 




h:m:s 










1 


03:33:17.9 


31:09:33.4 


1 


Bolo80 




2 


03:33:21.4 


31:07:30.7 


2,3 


B0I08I 




4 


03:33:16.3 


31:06:54.3 


4 






5 


03:33:01.9 


31:04:23.2 


5,6,8,9,11 


Bolo74 




7 


03:33:16.5 


31:07:51.3 


7 




03331667+3107548 


10 


03:33:27.3 


31:07:10.1 


10 






12 


03:43:56.5 


32:00:49.9 


12 


Bolol03 




13 


03:43:56.9 


32:03:04.8 


13 


Bolol04 




14 


03:44:43.9 


32:01:32.0 


14 


B0I0II6 


03444389+3201373, 












03444330+3201315 


15 


03:43:50.8 


32:03:24.2 


15 


Bolol02 


03435056+3203180 


16 


03:44:01.0 


32:01:54.8 


16 


B0I0IO6 




17 


03:43:57.9 


32:04:01.5 


17 


Bolol05 




18 


03:44:03.0 


32:02:24.3 


18 


Bolol07 




19 


03:44:36.8 


31:58:48.9 


19 


Bololl5 




20 


03:44:05.5 


32:01:56.7 


20,22 


BolollO 




21 


03:44:02.3 


32:02:48.5 


21 


Bolol07 




23 


03:43:37.8 


32:03:06.2 


23 


Bolo99 




24 


03:43:42.3 


32:03:23.2 


24 






25 


03:44:48.5 


32:00:30.8 


25 


Bololl7 




26 


03:43:44.4 


32:02:55.7 


26 


BololOO 




27 


03:25:35.9 


30:45:30.0 


27,34 


B0I08 


03253643+3045252 


28 


03:25:36.4 


30:45:15.0 


28,34 


B0I08 


03253652+3045070, 












03253643+3045252 


29 


03:25:38.8 


30:44:03.6 


29 


BololO 




30 


03:25:22.4 


30:45:10.7 


30 


Bolo5 




31 


03:25:25.9 


30:45:02.7 


31 


— 




32 


03:25:49.0 


30:42:24.6 


32,33 


Bolol3 




35 


03:27:39.1 


30:13:00.6 


35 


Bolo22 




36 


03:27:42.9 


30:12:28.5 


36 


Bolo23 




37 


03:27:48.4 


30:12:08.8 


37,38 


Bolo24 


03274767+3012043 


39 


03:27:38.1 


30:13:57.3 


39 


Bolo21 


03273825+3013585 


40 


03:27:39.9 


30:12:09.8 


40 






41 


03:29:10.4 


31:13:30.0 


41 


Bolo48 




42 


03:29:12.0 


31:13:10.0 


42 


Bolo48 




43 


03:29:03.2 


31:15:59.0 


43 


Bolo43 


03290375+3116039 


44 


03:28:55.3 


31:14:36.4 


44 


Bolo38 




45 


03:29:01.4 


31:20:28.6 


45 


Bolo42 


03290116+3120244, 












03290149+3120208, 












03290051+3120284 


46 


03:29:11.0 


31:18:27.4 


46 


Bolo49 




47 


03:28:59.7 


31:21:34.2 


47 


Bolo40 




48 


03:29:13.6 


31:13:55.0 


48 






49 


03:28:36.7 


31:13:29.6 


49 


Bolo29 


03283609+3113346 


50 


03:29:06.5 


31:15:38.6 


50 




03290642+3115348 


51 


03:29:08.8 


31:15:18.1 


51 


Bolo46 




52 


03:29:03.7 


31:14:53.1 


52 






53 


03:29:04.5 


31:20:59.1 


53 






54 


03:29:10.7 


31:21:45.3 


54 




03291063+3121469 


55 


03:28:40.4 


31:17:51.3 


55 


Bolo31 




56 


03:29:07.7 


31:21:56.8 


56 




03290773+3121575, 












03290704+3121578 


57 


03:29:18.2 


31:25:10.8 


57 


Bolo53 





Other associated source names 



bl-c (Hiranoet alj I19991J> 
[20051) 

bl-bS; bl-bN jHirano et alj 
JWalawender et alJ2005l) 
bl-d JHirano et alj H999bT) . Bl 
120051) 



Bl SMM2 (Walawenderetal. - 
I1999H) . Bl SMM1 
SMM3 (Walawenderetal. 



IRAS 03301 +3057, Bl SMM6 ( Walawen deret alj2005l) 
Bl SMM11 t Wala wenderet alj2005t) 

HH211 jMcCaughrean et alJll994tlChandler & Richerll2000l) : 
HH211 VLA 2jAvilaetalj2001l) 
IC348 MMS jEisloffel et aU2003l) 



none 



none 
none 



CXOPZ 15 iPreibisch & Zinne cker 2001) 



none 
none 
none 
none 
none 
none 



L1448 NW (Barsonv etal.1 119981; 


Bachiller & Cernicharol 


1986:Curiel etal.ll99C 


: Loonev et al. 


200bl) 


L1448 N A/B 


Barsonv et alj 119981 ICurieletalJ 



1 990 : lLoonev et alj 
jBadnlkr^^Cemichard 
jAnglada & Rodriguedl2002 



LI 448 IRS 
IRS2 



3A/3B 
VLA6 



L1448C JBarsonv et alj ~ l 19981 : iBachiller et al] 119911 : 
ICurielet all 19901) 

IRAS 03222+3034 L1448 IRS2 JO'Linger et alj [1999 
L1448 IRS2 VLA4 jAnglada & Rodriguez 2002h 



IRA S 03245+3002, L1 455 FIR4 jTapia et al.lfl997h . RN015 
FIR jRengel et aU2002I) 

none 

L1455 PP9 jTapiaetalll997l).G N 03.24.7 jMagakianll2003l) 
L1455 FIR1/2 ffapia et all 19971) 

none 

NGC1333 IRAS 4 A Jjennings et al] 119871) VLA25 
(Rodriguez et all 19971) . SK4 jSandell & Knedl2001h 
NGC1333 IRAS 4B Jjennings et aTl 119871) . SK3 
JSandell & Knee]|200j|j 

NGC 1333 SVS13 JStrom et alj [19761) . ASR 2 jAspin et al. 
1 19941) IRAS 03 259+3 105,HH7-1 1 M MS 1 JChini et al 
Il997ah . VLA 4a/b jRodriguez e7al] | 1997l) 
NGC1333 IRAS 2A Jjennings et all 19871) 
NGC1333 ASR 114 JAspin et al.1 119941). VLA 42/43 



(Rodriguez et all 19971) . SK24 JSandell & Knee 2001) 



NGC1333 ASR 32/33 jAspin et alj Il994h. IRAS 7 SMI/2 
JSandell & Knedl200ll).VLA 2 JSnell & Ballvll 19861) 

NGC1333 SK3 USandell & Knedl2001h 

NGC 1333 IRAS 4C Jjennings et alJ1987l:lLoonev et alj2000l) . 
SK5 JSandell & Kneel2001l). VLA29 iRodrfguez et alj|1999n 
NGC1333 SK6 JSandell & Kneel2001n, IRAS 03255+ 3103 
NGC1333 HH 7- 1 1 MMS 4 JChini et alj Il997bl) ASR6 
Aspinetal] 119941). 



vonHippel etal.ll988l) 



SK15 



2001), HH8 



NGC1333 SK16 JSandell & Knedl200l[) 
NGC1333 SK14 JSandell & Kned I200l1) . HH7-11 MMS 6 
JChini et all 1997bf) 

NGC1333 SK26 jSandell & Kjiee|[200l[) 
NGC1333 SK28 JSandell & Kn^l200l1) 

none 

NGC1333 SK29 JSandell & Knei200ll) 

NGC1333 SK33 JSandell & Knell200~R) 
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Table 5. continued. 



Source 


RA(J2000) 
h:m:s 


Dec(J2000) 
d:':" 


Paper I 


Bolocam 


2MASS 


Other associated source names 


58 


03:29:24.0 


31 


33:20.8 


58 


Bolo57 




none 


59 


03:29:16.5 


31 


12:34.6 


59 


Bolo51 




none 


60 


03:28:39.4 


31 


18:27.1 


60,73 


— 




none 


61 


03:29:17.3 


31 


27:49.6 


61 


Bolo52 




none 


62 


03:29:07.1 


31 


17:23.7 


62 


Bolo45 




NGC1333 SK18 jSandell & Knee||2001h 


63 


03:29:18.8 


31 


23:16.9 


63 


Bolo54 




NGC1333 SK32 jSandell & Kj^|20o!I) 


64 


03:29:25.5 


31 


28:18.1 


64 


Bolo58 




NGC1333 Per 4A3/4D (Ladd et al J 19*94) 


65 


03:29:00.4 


31 


12:01.5 


65 


Bolo41 




NGC1333 IRAS 4B1 (Jennings et alj Il987h. SKI 
fSandell & Knee 2~001h 


66 


03:29:05.3 


31 


22:11.3 


66 


— 




NGC1333 SK30 (Sandell & Knedl2001h 


67 


03:29:19.7 


31 


23:56.0 


67 


— 


03292003+3124076 


none 


68 


03:28:56.2 


31 


19:12.5 


68 


Bolo39 




NGC1333 MB0146 (Wilkins et alJ2004t) 


69 


03:28:34.4 


31 


06:59.2 


69 


Bolo28 




none 


70 


03:29:15.3 


31 


20:31.2 


70 


Bolo50 




NGC1333 SK22 (Sandell & Knee 2001) 


71 


03:28:38.7 


31 


05:57.1 


71 


Bolo30 




HH 340B (von Hirroel et all 19880 


72 


03:29:19.1 


31 


11:38.1 


72 


Bolo55 




none 


74 


03:28:32.5 


31 


11:07.7 


74 


Bolo25 


03283258+3111040 


NGC1333 38 (Lada et alJ1996T) 


75 


03:28:42.6 


31 


06:10.0 


75 


Bolo33 




none 


76 


03:32:17.8 


30 


49:46.3 


76 


B0I066 




IRAS 03292+3039, Bl SMM5 (Walawender et all2005h 


77 


03:31:21.0 


30 


45:27.8 


77 


Bolo65 




IRAS 03282+3035, Bl SMM7 (Walawender et alj20*05T) 


78 


03:47:41.6 


32 


51:44.0 
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Fig. A.l. The density structure of the inner region of the Class source used in the radiative-transfer modelling. 

Appendix A: Class models 

Below we provide details of the Class structure used in the radiative transfer modelling. 

A.1. Density structure 

We describe the geometry of the Class sources using a similar density structure to that of Whitnev et al. 1 2003). In summar y, the geometry consists of a central 
source of mass M,, radius R„, and temperature T c g = 4000 K surrounded by a canonical flared disc Isha kura & Sunyaev 1 9731) : 



p (f -, z )=po(^fexp(-I(£) 2 ) (A.1) 

where r is the radial distance in the midplane and z is the distance perpendicular to the midplane. The flare is imposed by a power-law dependence of the scale-height 
honr.h = hoir/R,)^. The disc mass is then used to define po via 

M di5C Q3 -a + 2) 

Po = -„ (A. 2) 

For all the models de scribed below we adopt the following disc parameters: /J = 1.25, a = 2.25, h = 0.01S,, = 0.01 M G , Rinner = 7.5 R Q , and R outer = 10 AU 
IWhitnev et al.120031) . 

The disc /star system is enshrouded by an infalling envelope of mass M env , which has a density structure described by a rotationally flattened infalling envelope 
4Ulrichll976D : 



An \ r 3 / \ po 



— + —5— (A.3) 
(Po r 



where M cn v is the mass infall rate, R c is the centrifugal radius, fi = cos 8 = zl V (z 2 + r 2 ) and /iq is the cosine of the polar angle of a streamline of infalling material 
as r — > oo, and is found by solving: 

fil + fioir/Rc - 1) - n(r/R c ) = (A.4) 

We fix the inner radius of the envelope at R mReI and set the outer radius to 5000 AU. We include narrow, fixed opening angle (5°) opening angle, low density (10~ 24 
g cm -3 ) cavities in the polar directions. The density structure of the inner region of this Class model can be seen in Figure IaTI 
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Fig. A.2. Model SEDs for Class objects. The flux of the object, in AF,\ (Wm 3 ) is plotted against wavelength in microns. 
Table A.l. Model parameters. The envelope mass (M env ), stellar mass (M„), and stellar luminosity (L*) are listed in solar units. 



M cnv (M ) 


M, (Mo) 


L, (L Q ) 


10.00 


3.33 


75.9 


3.00 


1.00 


14.9 


1.00 


0.33 


3.4 


0.50 


0.17 


1.3 


0.20 


0.07 


0.38 


0.10 


0.03 


0.15 


0.05 


0.02 


0.06 



A.2. Dust properties 

We adopt a standard MRN ( Mathis et al. 1977) grain-size distribution 

n(a) da cc a da (A. 5) 

where n{a) is the number density of grains with diameter a. We use minimum and maximum grain sizes of 0.005 /im and 1 fim re spectively. We assume th at the 
dust is composed of 62.5% 'astronomical' silicate and 37.5% graphite grains, with the usual 2:1 ratio of ortho- and para- g raphite (braine & Malhotra 1993]). The 
cross-sections and Mie phase matrices were computed directly from the refractive indices tabulated by Draine & Lee 1 1984J). A dustigas ratio of 1:100 was used for 
all models. 



